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ABSTRACT 

Context. Unevolved metal-poor stars constitute a fossil record of the Early Galaxy, and can provide invaluable information on the 
properties of the first generations of stars. Binary systems also provide direct information on the stellar masses of their member stars. 
Aims. The purpose of this investigation is a detailed abundance study of the double-lined spectroscopic binary CS 22876-032, which 
comprises the two most metal-poor dwarfs known. 

Methods. We have used high-resolution, high-S/N ratio spectra from the UVES spectrograph at the ESO VLT telescope. Long-term 
radial- velocity measurements and broad-band photometry allow us to determine improved orbital elements and stellar parameters for 
both components. We use OSMARCS ID models and the turbospectrum spectral synthesis code to determine the abundances of 
Li, O, Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Fe, Co and Ni. We also use the C0 5 BOLD model atmosphere code to compute the 3D 
abundance corrections, especially for Li and O. 

Results. We find a metallicity of [Fe/H] 3.6 for both stars, using ID models with 3D corrections of 0.1 dex from averaged 

3D models. We determine the oxygen abundance from the near-UV OH bands; the 3D corrections are large, -1 and -1.5 dex for the 
secondary and primary respectively, and yield [O/Fe] ~ 0.8, close to the high-quality results obtained from the [OI] 630 nm line in 
metal-poor giants. Other [o'/Fe] ratios are consistent with those measured in other dwarfs and giants with similar [Fe/H], although Ca 
and Si are somewhat low ([X/Fe]Sj 0). Other element ratios follow those of other halo stars. The Li abundance of the primary star is 
consistent with the Spite plateau, but the secondary shows a lower abundance; 3D corrections are small. 

Conclusions. The Li abundance in the primary star supports the extension of the Spite Plateau value at the lowest metallicities, without 
any decrease. The low abundance in the secondary star could be explained by endogenic Li depletion, due to its cooler temperature. 
If this is not the case, another, yet unknown mechanism may be causing increased scatter in A(Li) at the lowest metallicities. 

Key words, nuclear reactions, nucleosynthesis, abundances - Galaxy:halo - Galaxy:abundances - cosmology:observations - stars: 
Population II 
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* Based on observations made with the ESO Very Large Telescope at 
Paranal Observatory, Chile (Large Programme "First Stars", ID 165.N- 
0276(A); P.I. R. Cayrel). 



1. Introduction 

Extremely metal-poor (EMP) stars formed with the chemical 
composition of the gas in the early Galaxy, and constitute a 
unique source of information on the first generations of stars. 
Among EMP stars, a special place is held by the dwarfs, which 
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are not subject to the mixing episodes experienced by giants, 
thus enhancing their value as cosmological probes. 

In fact, among these stars the Li abundance appears to 
be constant whatever the stellar temperature or metallicity 
( Spit e & Spite 1982a b), the Spite plateau. The simplest inter- 
pretation of the plateau is that it represents the primordial Li 
abundance, i.e., it reflects the amount of Li formed in the first 
minutes of the existence of the Universe. If so, Li can be used 
as a "baryometer", a tool to measure the baryonic density of 
the Universe, since this is the only cosmological parameter upon 
which the primordial Li abundance depends upon. 

The independent determination of the baryonic density from 
the fluctuations of the Cosmic Microwave Background by the 
WMAP satellite dSpergel et al. II2003LI2007I) and other CMB ex- 
periments m easuring fluctuations on smaller angular sc ales, such 
as the VSA dRebolo et al. 1120041: iGrainge et al. 112003b . ACBAR 
dKuoetal.l 12004 and CBI dPearson et al. 1120031) experiments, 
implies a primordial Li abundance which is at least a factor of 
3-4 larger than that observed on the Spite plateau, creating a 
conflict with the traditional interp retation of the plate au. 

In Paper VII in this series dBonifacio et al. 1120071) we have 
investigated the Spite plateau at the lowest metallicities (down 
to [Fe/H]=-3.3) and found marginal evidence that at these low 
metallicities there could be an increased scatter or even a sharp 
drop in the Li abundance. It is therefore of great interest to ex- 
plore the Li abundance in stars of even lower metallicity. 

The star CS 22876-032 was identified in the fi rst paper re- 
portin g results of the HK objective-prism survey bv lBeers et al. I 
(119851) . who noted that it had the weakest Ca II K line in the 
low-metallicity sample, suggesting that it could be as metal- 
deficient as the record ho lder at that time, the giant CD -38°245 
JBessell & Norris 111984 . CS 22876-03 2 had already been ob- 
served in the objective-prism survey of ISlettebak & Brundage I 
(fl97l . who classified it as an A-type peculiar star and noted its 
weak and diffuse Balmer lines. Having assigned to this star a 
much earlier spectral type, they did not conclude that the weak- 
ness of the Ca II K line was indeed due to an extremely low 
metallicity. 

At the confere nce "Chemical and Dy namical Evolution of 
Galaxies" in 1989 dBonifacio et al. Ill 990h . P. Molaro announced 
that high-resolution spectra from the CASPEC spectrograph at 
the ESO 3.6 m telescope indica ted \Fe/U] ~ - 4.3 for CS 22876- 
032. However, just afterwards iNissen 1 dl989l) discovered, from 
simlar- resolution spectra, that the star is a double-lined spectro- 
scopic binary. The spectra acquired by Molaro were obtained at 
a singl e-lined phase, and the abun dance analysis of CS 22876- 
032 by Molaro & Castelli ( 1990) assumed that it was a single 
star. Thus, veiling was neglected, the adopted temperature was 
too low, and the measured [Fe/H] was therefore a lower limit to 
the metallicity of the system. 

Although CS 22876-032 is relatively bright (V= 12.84) for 
an EMP star, it took another ten years before a sufficient num- 
ber of high-resolution spectra had been accumulated to allow 
to determine of the orbital parameters of this system, and to 
perfor m a consistent chemical analysis. iNorris. Beers. & Ryan I 
( 2000) found the orbital period to be 424.7 days and the metallic- 
ity of the system [Fe/H] =-3.71. In spite of the upward revision 
of the metallicity, partly due to the different solar Fe abundance 
assum ed (log e(X) = 7.50 instead of 7.63 in lMolaro & Ca stelli 
I1990I) . the two stars in CS 22876—032 remain the most metal- 
poor dwarfs known. 

Thus, the CS 22876-032 system constitutes a unique fos- 
sil, recording the chemical composition of the early Galaxy. 
Moreover, it allows a measurement of the Li abundance which 



probes the Spite plateau at a lower metallicity than any other 
known dwarfs. Note that, despite its the extremel y low iron 
abund ance ([Fe/H] =-5 .4), the star HE 1327-2326 dFrebel et al. I 
2005) has very high C, N and O abundances, so its global metal- 
licity, Z, is considerably higher than that of CS 22876-032. It has 
also been shown recently that this star is most likely a slightly 
evolved subgiant, not a dwarf. 

In this paper we use high-resolution, high-S/N ratio spectra 
from the ESO Kueyen 8.2m telescope and the UVES spectro- 
graph to improve the orbital solution and perform a complete 
chemical analysis o f the two stars that comprise C S 22876-032. 
With respect to the Norris, B eers. & Ryan I d2000b analysis, our 
superior S/N ratio and larger spectral coverage permit measure- 
ment of abundances for many mo re elements, and, most i mpor- 
tantly, for both components; the INorris. Beers. & Ryan I ([2000) 
analysis of the secondary star was limited to Fe. 
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Fig. 1. Upper panel: Radial velocities of CS 22876-032 (filled 
circles: A; open circles: B). Open diamonds: Data from Norris 
et al. (2000). The curves show the orbital solution (P = 425 d, 
e = 0.14) for Star A (solid) and B (dashed). Dot-dashed hori- 
zontal line: Centre-of-mass velocity of the system. Lower panel: 
Residuals from the fit. 



2. Observations and data reduction 

Spectroscopic observations of the CS 22876-032 were carried 
out with the UV- Visual Echelle Spectrograph (UVES, Dekker 
et al. 2000) at the European Southern Observatory (ESO), 
Observatorio Cerro Paranal, using the 8.2 m VLT-Kuyen tele- 
scope on 2000 July 19, 20, August 3, 11, and October 17, 20, 
and 2001 November 7, 8 and 9, covering the spectral region 
from 300.0 nm to 1040.0 nm. Most of the observations were 
made with a projected slit width of 1" at a resolving power 
A/6A ~ 43 000. The spectra were reduced in a standard manner 
using the UVES reduction package within the MIDAS environ- 
ment. The signal-to-noise ratio per pixel varies from 25 at 312.0 
nm, 50 at 330.0 nm up to 150 or higher above 410.0 nm. 
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Table 1. Radial-velocity observations of CS 22876-032. For each velocity, cr is the estimated error and ( O-C) the residual from the 
orbital fit. *: H a velocities; omitted from solution. 



Date 


HJD-2,400,000 


Phase 


V A 


o-V A 


(O - C) A 


V B 


<tVb 


(O - C) B 


Source 


(days) 




km s _1 


km s _I 


km s 


km s" 1 


km s 


km s _I 






i no ^ Ca« A 
lyoj oep o 


40.51 J. 1 


n AAn 
0.66y 


01 n 

-83. y 


n 
Z.O 


1AA 

-3.44 








wd r> nn 


i no £ r~w. 1 a 
IVoJ Dec 10 


A £.A 1 A H 

4o41u.U 


n om 

o.yu / 


OA O 


n 
Z.O 


1 -70 
1. /Z 








mtj r> nn 


iy88 oep 2/ 


A HA 1 1 H 
4/431. / 


n on-7 

o.2y / 


nQ c 

-yj.3 


1 n 
1 .0 


1 A A 

1.44 








Mcyo 


iysy oep i i 


4/ 


n 1 0/1 
U. 1Z4 


1 fY7 O 
-10/. 8 


1 n 
1 .1) 


n aq 
0.63 


-7/1 O 
- /4.8 


1 n 
1 .U 


1 A O 

1.48 


vtd d nn 


iysy uct is 


AHQ 1 A 
4/014.0 


n 1 nn 

0.1 yy 


1 n/i q 
-104.3 


1 n 
1.0 


n n/i 
-0.94 








mo 15 nn 
JNBKOO 


iysy uct io 


A HQ 1 C 1 
4/0l J. 1 


n on 1 
0.2U1 


1 01 1 
-103.1 


1 n 
1 .0 


n no 
0.0/ 


OIO 

-81.2 


1 n 
1.0 


n oa 
0.86 


Mil nn 
JNBKOO 


i non 1 n 
lyoy (Jet 1 / 


AHQ 1 


n on/i 
0.2U4 


1 no 1 
-1U2.1 


1 n 
1 .0 


n o"7 
0.6/ 








JNBKOO 


1 ooo rv.,. A 
lyoy Dec 


AHQ6.H (\ 


0.3ZZ 


-yi.i 


1 n 
1 .U 


1C 

-0.33 








mtj d nn 
JNBKUU 


i nnn c . . 1 

iyyo oep 1 


/I O 1 1 Q H 


Oder 


1 ni n 
-103. 


1 n 
1.0 


n oi 
0.83 








/ ' A i TIT ' / ' 

CAorJbC 


i nnn c a « o-7 
iyyo oep LI 


A O 1 AO 1 

4ol02. 1 


(i on 
0.01 / 


1 nn A 

-iuy.6 


1 n 
1 .0 


n n 1 

-o.y 1 


-7A O 
- /O.Z 


1 n 
1 At 


n on 
-O.ZU 


MT3 r> nn 
JNBKUU 


1 (If! 1 A /fur T1 

lyyi May Z3 


/i /inn n* 


n c-7-7 
0.3 / / 


0"7 C 

-8 /.3 


n 
2.0 


H 1 A 

- 1.34 








tMMl 


iyyi Aug 22 


AO A c\n c 

48490.5 


0.790 


-87.0 


2.0 


-0.47 








\ * 1 1 [ t T W T 

JNBKOO 


i nn 1 c . . ~ o i 
lyy 1 oep 21 


A oin O" 

48520. / 


n oa 1 
0.861 


n/i n 

-y4.o 


1 n 
1.0 


n ac 
-0.63 








tMMl 


i nno a ~ n 

iyy2 Aug 1 / 


A OOO 1 O 

46531.6 


n A/i n 
0.640 


-70 A 

-/8.6 


1 n 
1.0 


1/11 
1.41 


1 n-7 a 
-1U/.6 


1 n 
1.0 


n 1 
-0.1Z 


mo nn 


i nno r~r 1 o 

iyy2 Dec 12 


/lOnAO" C* 

4890 I.J 


on 
0.912 


nn n 

-yo.y 


n 
2.0 


O Q Q 

8.33 








eMMI 


i nnA a ...» o" 

iyyo Aug / 


30303.1 


n ooi 
0.03 / 


1 1 n 
-1 10.2 


1 n 
1.0 


n 00 
-0.Z8 


-7A O 
- /O.Z 


1 n 
1.0 


1 ?c 
-1.33 




i nn-7 a >.~ oq 
iyy / Aug Z3 


c/vol O 

50683.8 


n nc7 
0.93 / 


1 r\A 
-1U4.2 


1 n 
1 .0 


n c-7 
-0.3 / 


O 1 T 

-81.3 


1 n 
1.0 


n oc 
0.Z5 




i nno a ...» it 

iyyo Aug 12 


c 1 mo q 
31038. 3 


n 101 
0. /8 / 


00" 
-8 /.8 


1 n 
1.0 


1 An 
-1.4/ 


1 m 1 
-1U1.1 


1 n 
1.0 


n 

-0.36 


mtj nn 


i nnn t..i on 

iyyy Jul 2y 


31388. / 


n A i o 
0.612 


"7-7 n 
-/ /.0 


n 
2.0 


on 

z.sy 








Mi] nn 


1999 Sep 23 


51444.6 


0.744 


-84.3 


1.0 


-0.93 


-105.5 


1.0 


-1.71 


NBR00 


2000 Jul 19 


51744.8 


0.451 


-84.9 


0.5 


-0.52 


-102.9 


0.3 


-0.22 


VLT 


2000 Jul 20 


51745.8 


0.453 


-83.7 


0.4 


0.56 


-102.3 


0.8 


0.51 


VLT 


2000 Aug 3 


51759.8 


0.486 


-82.6 


0.3 


0.13 


-103.9 


0.5 


0.59 


VLT 


2000 Aug 1 1 


51767.9 


0.505 


-82.3 


0.3 


-0.31 


-105.3 


0.3 


0.01 


VLT 


2000 Oct 17 


51834.7 


0.662 


-80.1 


0.3 


0.22 


-107.0 


0.4 


0.14 


VLT 


2000 Oct 20 


51837.6 


0.669 


-79.9 


1.0 


0.56 


-105.6 


1.4 


1.38 


VLT 


2001 Nov 7 


52220.6 


0.571 


-80.6 


0.3 


-0.35 


-107.4 


0.3 


-0.18 


VLT 


2001 Nov 8 


52221.6 


0.573 


-80.2 


0.3 


0.01 


-107.2 


0.3 


0.06 


VLT 


2001 Nov 9 


52222.6 


0.575 


-79.4 


0.7 


0.78 


-106.6 


0.9 


0.69 


VLT 



Sources: NBR00: iNorris. Beers.&Rvanl d2000l. and references therein); MC90: iMolaro & Castelli I dl990h : CASPEC: Unpublished 
velocity from a Caspec spetrum; EMMI: Unpublished velocities from EMMI spectra; VLT: This paper. 



Table 2. Orbital elements of CS 22876-032 



Parameter 


This paper 


Norris et al. (2000) 


P (days) 


424.81 +0.37 


424.71 ±0.60 


T a (HJD-2,400,000) 


48579.8 +7.1 


48576.4 ± 13.5 


e 


0.143 ±0.013 


0.12 ±0.03 


w (deg) 


148.1 ±6.4 


144.96 ± 12.4 


Vo(km s- 1 ) 


-93.11 ±0.13 


-93.36 ± 0.28 


K A (kms- 1 ) 


15.04 ± 0.26 


15.13 ±0.51 


K B (kms" 1 ) 


16.51 ±0.27 


17.06 ±0.56 


M A sin 3 i(M e ) 


0.701 ± 0.021 


0.76 ± 0.04 


M B sin 3 i(M Q ) 


0.639 ± 0.019 


0.68 ± 0.04 


Mb/Ma 


0.911 ±0.022 


0.89 ± 0.04 


cr (km s~') 


1.01 


1.60 



3. Revised orbital parameters 

We derived radial velocities from the UVES spectra by fitting a 
Gaussian to several unblended spectral lines within the IRAlQ 
context. Table Q] contains the radial velocities and the 1-cr er- 

1 IRAF is distributed by National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Foundation. 



rors estimated from the dispersion of the measurements of dif- 
ferent stellar lines. We also list other velocity data given by 
INorris. Beers. & Ryan I (2000) and references therein, or which 
we have measured from previously unpublished CASPEC or 
EMMI spectra of this system. 

Here and in the rest of the paper, we denote the more massive 
and luminous primary star as Star A, the secondary as Star B. 

Our new radial-velocity measurements of CS 22876-032 ex- 
tend the time coverage of its orbit considerably and permit im- 
pr ovement of the orbital element s relative to those by published 
bv lNorris. Beers. & Ryan I (|2000). The computed orbital param- 
eters are given in Table |2j Fig. [TJ compares the observed radial 
velocities of both stars with the curves predicted from these or- 
bital elements. 

Note especially the improved mass ratio, Mb /Ma = 0.91 + 
0.02, which provides stringent constraints on the stellar param- 
eters for the two components of the binary as discussed below. 
Note also that the orbital eccentricity is the lowest found among 
halo spect roscopic binaries with periods in the rang e 100-2000 
days (see lLafham et al. ll2002l;lGoldberg et al. 112002b . While this 
might be a hint that tidal interaction has been strong in this sys- 
tem, perhaps in the pre-main-sequence phase, the separation of 
the stars has been so large throughout their main-sequence life 
that this is unlikely to be the cause of the apparent Li depletion 
we find in star B (see Sect. 16. It . 
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4. Model Atmospheres 

4. 1 . One-dimensional models 

Our analysis used OSMARCS ID LTE model atmospheres 
(iGustafsson et al. I Il975_ iPlez et al.l Il992t lEdvardsson et al. I 
Il993t lAsplundetal.l Il997t IGustafsson et al. I I2003I) and th e 
turbospectrum spectral synthesis code ([Alvarez & Plez II1998I) . 
Models were interpolated in pre-computed grids for a metal- 
licity of [Fe/H] = -3.6, since this was the final iron abun- 
dance (see § I5.2I ). and with an a-element enhancement of 
[q/Fel = +0.4 dex. W e adopted solar abundances from 
iGrevesse & SauvaTI d2000t) . with the exception of O, for which 
we adopted log e(0) - 8.72, based on 3D mo del atmospheres 
dLudwig & Steffenl I2007I: ICaffau et al. Il2007bl) . The code tur- 
bospectrum is used to determine ID element abundances in each 
component of the binary, either via spectrum synthesis or by 
comparing the observed equivalent widths of different stellar 
lines with the theoretical curves of growth (see § !5.2l i. 
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Fig. 2. Average temperature profile of 3D CO s BOLD atmo- 
spheric models (solid lines) compared to IDlhd (dashed) and 
ID MARCS models (dotted) in stars A (bottom panel) and B 
(top). 



4.2. Three-dimensional models 

In this work we also consider two 3D model atmospheres, which 

s r 1 

have been computed with the CO BOLD code dFrevtag et al. I 

l2002t IWedemever et al. I [2003), one for each star. The atmo- 
spheric parameters are close to those observed for the two stars: 
T e « Aog g /[Fe/H]: 6550/4.50/-3.0 (A) and 5920/4.50/-3.0 (B). 
Each model consists of a representative set of snapshots sam- 
pling the temporal evolution of the photospheric flow at equal 
intervals in time. The total time intervals were 2400 s for the 
warmer star A, and 9500 s for the cooler star B. These time inter- 
vals should be compared to the convective turn-over timescales. 
From the hydrodynamical point of view, typical timescales in the 
models for both components are not much different from that in 
a solar model, where the convective turn-over timescale amounts 
to about 500 s. Thus, we sample about five turn-over timescales 
for the hotter and almost twenty for the cooler component. 

The comparison of 3D vs. ID models depends on which par- 
ticular ID model is chosen. We compared each of our 3D models 
(hereafter denoted as (3D), obtained from the mean temperature 



and pressure structure of the full 3D model), to a corresponding 
standard hydrostatic ID model atmosphere (hereafter denoted 
as IDlhd)- The (3D) model is a temporal and horizontal average 
of the 3D structure over surfaces of equal (Rosseland) optical 
depth. It is only dependent on the particular way the 3D model 
is averaged. 

The IDlhd model is calculated with a ID atmosphere code 
called LHD. It assumes plane-parallel geometry and employs the 
same micro-physics (equation-of-state, opacities) as CO s BOLD. 
Convection is described by mixing-length theory. Somewhat ar- 
bitrary choices to be made relate to the value of the mixing- 
length parameter, which formulation of mixing-length theory 
to use, and ho w turbulent pressur e is treated in the momentum 
equation; see ICaffau et al.l ([2007) for further details. As usual, 
in the spectral synthesis of the ID models a value of the micro- 
turbulence has to be adopted. For ID as well as 3D models the 
spectral synthesis calculations were performed with the spec- 
trum synthesis code Linfor3D0. 

4.3. 3D corrections 

There are two main effects that distinguish 3D from ID models, 
the average temperature profile and the horizontal temperature 
fluctuations. We quantify the contribution of both effects by in- 
troducing the 3D correction as: 3D - IDlhd. 

The average temperature profile provided by a hydrodynam- 
ical simulation is different from that of a ID atmosphere assum- 
ing radiative equilibrium. This effect is shown in Figf2j where 
the 3D average temperature profile, plotted as a function of the 
pressure, is compared to the profiles from IDlhd and MARCS 
models. As evident from the plot, the (3D) temperature profile 
is cooler than both ID models in the outer photospheric layers 
for both of th e stars. This often-enc ountered effect in metal-poor 
atmospheres dAsplund et al. I I 1999) is particularly important for 
the oxygen abundances derived from OH molecules, as the dif- 
ference is largest precisely in the layers where these lines are 
formed. The result is that the oxygen abundances become lower 
in the 3D formulation than in the ID; we quantify this effect 
through the 3D correction as defined above. 

5. Chemical analysis 

5.1. Stellar parameters 

The atmospheric parameters of each star in the CS 22876- 
03 2 system wer e estimated f rom the photometric data available 
in INorris et aO (Il993l) and iPreston et aTl (Il99lh . from whom 
we adopt V = 12.84, B - V = 0.397, U - B = -0.255 
with uncertainties of 0.01, 0.02 and 0.01 respectively. W e adop t 
E(B - V) = 0.00 ± 0.01 from INorris. Beers. & Ryanl (2000) 
and ISchuster et aTl dl993l I19961) . We also extracted, from the 
2M AS 50 database, K = 11 503 ±0.035 and J = 11. 800 ±0.041. 
The equations derived by Carpenterl d200ll) to transform from 
2MASS magnitudes to the homogenised photometric system of 
iBessell & Brett I dl988l) were then applied. 

From the above information it is possible to estimate the 
reddening-corrected colours (U - B)o, (B - V)o, (V - K)q and 

2 more information on Linfor3D can be found in the following link: 
http://www.aip.de/~mst/Linfor3D/linfor_3D_manual.pdf 

3 The Two Micron All Sky Survey is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 
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Fig. 3. The two components of CS 22876-032 on a 14-Gyr 
isochrone by Chieffi & Limongi (priv. comm.) for Z = 1CT 6 
(solid line). For comparison, isochrones for 14 Gyr and Z = 1CT 5 
(dashed) and 12 Gyr and Z = 1CT 6 (dotted) are also shown. 



(J - K)q, which we use to derive the parameters of both compo- 
nents of the binary by comparing with theoretical isochrones. We 
have chosen the isochrone of Chieffi & Limongi (private com- 
munication) for 14 Gyr and metallicity Z = 10~ 6 , from which 
one can compute composite colours from pairs of two models 
that lie on that isochrone. Thus, the stellar parameters can be de- 
rived from the best fit to the observed colours that also satisfy 
the mass ratio determined from the orbital solution. 

We note that the isochrones of Chieffi & Limongi use 
the colour transformations based on ATLAS model atmo- 
spheres, for this temperature r ange, and the synthetic colours of 
iBessell. Castelli & Plez I (fl998h . 

The result is shown in Fig.[3j which corresponds to a primary 
star with 7^ ,a = 6500 K and log g& = 4.4 and a secondary star 
with r effiB = 5900 K and log g B = 4.6. We checked that a 
change of AAge = -2 Gyr translates into a change of +25 K and 
-0.01 dex for the T e s t A,B and log gA,B respectively, whereas a 
variation of the metallicity of A log Z = + 1 dex does not have 
any impact on the resulting stellar parameters (see Fig. [3j. 

The uncertainties of the stellar parameters were estimated us- 
ing Monte Carlo techniques. We injected noise in the seven ob- 
served quantities, V, B-V, U-B, K, J, E(B-V), and M B /M A fol- 
lowing Gaussian distributions with standard deviations equal to 
the errors of these quantities. From these distributions we com- 
puted a set of five variables, {U-B) Q , (B-V\ h (V-K) , (J-K) , 
and Mb/Ma for the 10,000 samples. Then we found the best fit 
to each of these set of variables via a x 1 minimisation, defining 
X 2 = 2f=iC/i,obs - Amod) 2 , Aobs being the "observed" value for 
each Monte Carlo simulation and //. mo d the value extracted from 
two pairs of points in the theoretical isochrone. 

The results of these simulations for the effective temperature 
and surface gravity of both components are shown in Figs. [4] 
and [5] respectively. The lowest contour encloses roughly 95.4% 
of the 10,000 Monte Carlo events, analogous to 2cr for a normal 
distribution. From these simulations we adopted an error, at the 
2a level, of AT effA = lOOKand AT eB , B = 150 K for the effective 
temperature, and Alog g^,B = 0.1 dex for the surface gravity. 

For single stars, t he wing s of Ha is also a very good 
temperature indicator dCavrell fl988t iFuhrmann et aTl 11993; 
Ivan t Veer-Menneret & Megessier II 19961: iBarklem et al"ll2002l) . 
Adopting the broadening theory of IBarklem et al. I d2000l) . we 
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Fig. 4. Distribution of effective temperatures for CS 22876-032 
A and B obtained by Monte Carlo simulations, comparing the 
observed colours with those from the isochrone in Fig.[3]for the 
observed mass ratio. The lowest contour encloses 95.4% of the 
10,000 simulations. 



computed Ha profiles for several effective temperatures, using 
turbospectrum. Unfortunately, all our UVES spectra were ob- 
tained near maximum line separation, and the velocity difference 
is of the order of ~ 30 km s _1 , which precludes separation of the 
individual Ha profiles. Therefore, we had to compute a compos- 
ite spectrum in order to match the observed profile. 

Fig.|6]compares this synthetic Ha profile with the observed 
profile for several combinations of effective temperatures T e g ^ 
+ T e g t B- We did not vary T e g t B because the Ha absorption line of 
the cooler star B is weaker and severely veiled by the flux of star 
A, so the combined profile is not sensitive to changes in T e ff >B . 
This comparison seems to confirm our estimate of the effective 
temperature from the colours, suggesting as well that it is on the 
same scale as the Ha-based temperatures. 

5.2. Stellar elemental abundances 

Most of the elemental abundances were determined from equiv- 
alent width measurements of selected unblended lines. These 
were made wit h an automatic line-fi tting procedure based on the 
algorithms of Cha rbonneau I (1 19951) . which performs both line 
detection and gaussian fits to unbl ended fines. Th e imple men- 
tation is the same as described in Francoi s~et al. I (J2003). The 
equivalent widths were then corrected for the appropriate veil- 
ing factors and provided as input to turbospectrum to determine 
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Tabl e 3. Element abundances of CS 22876-032. The Solar O abundance is adopted fr om the 3D measurement of atomic lines of 
iLudwig & Steffenl d2007l) ; the other Solar abundances from lGrevesse & Sauvall d2000t) . [X/H] and [X/Fe] are LTE values; [X/Fe] 
refers to Fei for OH and neutral species, to Fen for ionised species. The 3D abundance corrections A3D-1D were determined with 
turbospectrum from the (3D) and ID models, except for O and Li, where a full 3D analysis was performed and we give two values: 
The abundance correction 3D - (3D) ("(3D)") and "ID" = 3D - IDlhd (see text), cr is the standard deviation of the results from 
the n lines (next column; if n — 1, the wavelength of the line in nm is given). See text for details on the measurement of Li, O, Sc, 
and Co. 



Species 


log e^Ajo 






A3D-1DA 


0"A 


"a 


[A/njg 


TX /Fein 


^3D-1D,B 




"B 


Li 1 




2.22 




(3D):-0.17 
1D:-0.19 


0.01 


670.8 


1.75 




(3D):-0.02 
lD:-0.29 


0.04 


670.8 


O(OH) 


8.72 


-1.52 


2.14 


(3D):-0.64 
1D:-1.49 


0.04 


4 


-1.75 


1.81 


<3D}:0.00 
lD:-0.92 


0.09 


9 


Nai 


6.33 


-3.69 


-0.03 


-0.01 


0.03 


2 


-3.79 


-0.22 


0.02 


0.05 


589.5 


Mgi 


7.58 


-3.11 


0.55 





0.18 


9 


-3.14 


0.43 


0.09 


0.24 


9 


All 


6.47 


-3.89 


-0.23 


-0.02 


0.22 


2 


-3.75 


-0.18 


0.04 


0.22 


2 


Sil 


7.55 


-3.75 


-0.09 


-0.01 


0.02 


390.5 


-3.48 


0.09 


0.07 


0.04 


390.5 


Cai 


6.36 


-3.65 


0.01 


-0.07 


0.02 


422.6 


-3.68 


-0.11 


-0.04 


0.04 


422.6 


Can 


6.36 


-3.68 


-0.18 


0.10 


0.14 


3 


-3.81 


-0.54 


0.23 


0.1 1 


2 


Sc n 


3.17 


-3.73 


-0.22 





0.18 


2 


-3.57 


-0.30 


0.04 


0.11 


2 


Tin 


5.02 


-3.25 


0.25 


-0.01 


0.34 


19 


-3.45 


-0.18 


0.09 


0.23 


12 


Cri 


5.67 


-3.79 


-0.13 


-0.07 


0.12 


5 


-3.92 


-0.35 


-0.08 


0.35 


4 


Crn 


5.67 


-3.38 


0.12 


0.08 


0.10 


313.2 


-3.30 


-0.03 


0.24 


0.30 


313.2 


Mnn 


5.39 


-4.01 


-0.51 


0.09 


0.17 


344.1 


-3.86 


-0.59 


0.14 


0.37 


344.1 


Fei 


7.50 


-3.66 





-0.12 


0.12 


38 


-3.57 





-0.07 


0.21 


37 


Fe n 


7.50 


-3.51 





0.09 


0.25 


13 


-3.27 





0.18 


0.31 


9 


Co 1 


4.92 


-2.91 


0.75 


-0.13 


0.07 


7 


-2.92 


0.65 


-0.07 


0.08 


4 


Nil 


6.25 


-3.49 


0.17 


-0.19 


0.18 


12 


-3.45 


0.12 


-0.08 


0.16 


12 



Table 4. Abundance errors in CS 22876-032. Aj dr , Ai og g, and A^ are the abundance changes caused by changes in T e s of 100K (A) 
or 150K (B), in \ogg by 0.1 dex, and by 0.5 km s _1 in the microturbulence velocity Other column headings and comments as in 
Table[3] 
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the abundances. The detailed line-by-line abundances, together 
with observed EWs, veiling factors, and atomic data for both 
components can be found in Table I A. fl 

The mean abundances for each element, listed in Table[3] are 
computed using the adopted stellar parameters derived in § 15.11 
and a microturbulence of ^ t = 1-5 kms . We have only per- 
formed full 3D computations with Linfor3D for the Li doublet 
and the OH lines. A full 3D analysis for the hundreds of lines 
involved in this work is a considerable computational task, well 
beyond the scope of the present paper. However, we used the 
(3D) models as input to turbospectrum to estimate the expected 
corrections due to the different average temperature profiles of 
the 3D models. From the full 3D spectrum synthesis performed 
for Li and OH, we expect this to be close to the true 3D correc- 
tion for star B, while we expect significant contributions from 



the temperature fluctuations in star A. In the following we re- 
fer to ((3D) - IDmarcs) to represent the difference between the 
abundance found by turbospectrum using a (3D) model and that 
found using a MARCS ID model. This is to distinguish them 
from the true 3D corrections. 

5.2.1. Veiling corrections 

In a double-lined spectrum the strength of each spectral line, in 
particular those of the fainter component, is reduced by veiling 
from the continuum flux of the other star. Thus, the measured 
equivalent widths for each spectral component must be corrected 
for this veling effect, in order to obtain the intrinsic values. The 
corrected equivalent widths can be estimated by multiplying the 
observed EWs by veiling factors, faj, which solve the equation 
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Fig. 6. Computed H a profiles for three sets of effective tempera- 
tures, 6300+5900, 6500+5900 and 6700+5900, and normalised 
to the level of the observed UVES spectrum at 645.0 nm. 



1/Aa + = 10 and where fa B /fu is the primary-to- 

secondary luminosity ratio. The values f A j are wavelength de- 
pendent and can be estimated theoretically by computing the 
flux of each stellar component, taking into account the ratio of 
the stellar radii. Thus, the luminosity ratio can be expressed as 
Ab//aa = F AA /F A , B x (R A /R B ) 2 , where F AJ and R, are the flux 
and radius of each star. 



For consistency wi th the isochrone colour s, we use version 
9 of t he ATLAS code dKuruczlll993a[. l2005al) in its Linux ver- 
sion (Sbo rdone et alj|2004t ISbordond feOOS) to compute model 
atmospheres and fluxes for each star by adopting the r e ff and 
log g derived in § 15.11 We used the "N EW" Opacity Distribution 
Functions (Castelli & Kuruczl 120031) . with 1 km s micro- 
turbulence, a mixing-length parameter ffMu of 1.25, and no 
overshooting. The formulation of the mixing length is different 
between MARCS and ATLAS; that used in our MARCS models 
roughly corresponds to <tmlt ~ 1.1 in the ATLAS formulation. 

The veiling factors thus computed differ from the 
MARCS +TURBOSPECTRUM ones on average by 0.02 which, for a 
line on the linear part of the curve of growth, translates into a 
difference in abundance of 0.009 dex, which is negligible within 
the accuracy of our analysis. This allowed us to compute quickly 
veiling factors for different combinations of the parameters of 
the two stars, which were then used to estimate the associated 
uncertainties. For consistency, we used the ATLAS veiling fac- 
tors throughout. 

The ratio of the stellar radii was extracted from the theoret- 
ical isochrone, being Ra/Rb = 1.4. The derived veiling factors 
lie in the range f AA = 1.29 - 1.37 (f A , B = 4.47 - 3.70) in the 
spectral region AA3Q0. 0-700. nm. These estimates also com- 
pare well with those used by Nor ris. Beers. & Ryan I d2000l) . al- 
though they adopted fixed values for large spectral regions. In 
particular, they used f AA = 1 -28 and f AtB = 4.60 for all Fe lines 
between 370 and 440 nm, whereas we used f AA = 1.28-1.31 
and f A B — 4.19 - 4.58 in that spectral region (see Table lA.ll . 

In addition to the ID veiling factors computed for the ef- 
fective temperatures of both stars, we also calculated the (3D) 
veiling factors using the (3D) atmospheric models (whose tem- 
perature structure is different from that of the OSMARCS ID 
models). These veiling factors were adopted to properly correct 
the EWs given as input to the turbospectrum code to estimate 
the ((3D) - IDmarcs) abundance corrections (see § I5.21 i. In ad- 
dition, we note that the 3D - IDlhd correction (see § |4.2t , com- 
puted only for Li and OH lines, does not consider different veil- 
ing factors for ID and 3D models. In this case, we computed 
only the 3D veiling factors using the continuum flux provided 
by the full 3D model of each star in the spectral region close to 
these lines. Thus, these 3D veiling factors were applied to the 
observed EWs, and the resulting EWs were used to compute the 
3D - IDlhd corrections reported in Table [3] 

5.2.2. Uncertainties in the abundance analysis 

The abundance measurements are dependent on the model pa- 
rameters, i.e. effective temperature, surface gravity and micro- 
turbulence. However, in the analysis of a spectroscopic binary, it 
is not possible to avoid the influence of the veiling factors on the 
error estimates. The veiling factors depend on the effective tem- 
peratures and surface gravities of both stars. Therefore, in order 
to estimate the sensitivity of an element's abundance to a given 
stellar parameter, one should also estimate how veiling factors 
change when one of the stellar parameters of each star is mod- 
ified. Thus, we have also computed ID veiling factors for four 
pairs of models T e ff A flog gA, T e ff B /\og g B according to the er- 
rors of the stellar parameters (see § 15. IK by changing one stellar 
parameter and fixing the three remaining parameters. The un- 
certainties on the elemental abundances due to the errors of the 
different model parameters are listed in Table|4] The uncertainty 
of the microturbulence was assumed to be 0.5 km s _I . 

The errors computed from the dispersion of the line mea- 
surements and the signal-to-noise ratios are listed in Table [3] 
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Table 5. Li abundances in CS 22876-032. A(Li)iD,NLTE,DC includes corrections for depletion by diffusion (see text). 



Component 


EW obs (pm) 
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A(Li) 1D 


A(Li)lD,NLTE 
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A3D-<3D> 


A3D-1D 
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1.32 


1.36 


2.22 


2.18 


2.18 


-0.17 


-0.19 


StarB 


0.42 


3.74 


1.75 


1.77 


1.84 


-0.02 


-0.29 



We used the Cayrel formula (ICavrellll988l) to estimate the er- 
rors of the observed EWs. Due to the high S/N of the spectra 
we have obtained, these errors are typically ^0.1 pm in most 
of the spectral region covered, except for the blue spectrum at 
310.0-320.0 nm where the S/N ratio drops significantly. Thus, 
the dispersion of the measurements for elements with lines in 
this range is larger than for lines above 400 nm. In Table [3] we 
give the larger of these two estimates for each element in each 
star. For Li, the error due to the uncertainty in the continuum 
level has been computed from a Monte Carlo simulation, by in- 
jecting noise corresponding to the actual S/N ratio near the Li 
line in the best-fit synthetic spectrum. In each case, the S/N ra- 
tio was estimated taking into account the corresponding veiling 
factors. 

5.2.3. The iron abundance 

We made a careful selection of 38 reliable Fe I lines in star A 
and 37 in B, taking into account the radial velocity separation of 
the two sets of lines in the double-lined spectrum. It is reassur- 
ing that we find the same [Fe/H] for both stars, within the errors. 
Note t hat in our analysis, in contrast with lNorris. Beers. & Rvanl 
(2000), this was not imposed a priori, and thus supports our de- 
termination of atmos pheric parameters. Our value of [Fe/H] is 
very close to that by iNorris. Beers. & Ryan I (|2000|) despite the 
very different effective temperatures adopted. This results from 
the fact that different effective temperatures also imply different 
veiling factors, which must be factored in. 

Reviewing the ((3D) - IDmarcs) corrections, we see that 
[Fe/H] is reduced in 3D, and the difference between the two stars 
increases. We expect the true 3D corrections of the primary star 
to be largely due to horizontal temperature fluctuations, there- 
fore probably larger than our ((3D) - IDmarcs) correction, and 
we predict that a full 3D-LTE synthesis will not improve the 
agreement of [Fe/H] between the two stars. Our best estimate 
of [Fe/H] of the system is still the ((3D) - IDmarcs) corrected 
value for star A: [Fe/H]=-3.78, which might be further reduced 
by a full 3D-LTE synthesis. This confirms that the stars in CS 
22876-032 are the most metal-poor dwarfs known to date. 

We note that ionisation equilibrium is not achieved in either 
star. In both, the abundance derived from the Fe n lines is larger, 
for star B by a factor of two. We note, however, that the Fe n 
abundance shows a very large scatter in both A and B (0.25 and 
0.31 dex, respectively) so that, within errors, the Fei and Fen 
abundances remain compatible. The number of Fe 11 lines mea- 
sured is very large for stars of this metallicity (13 in star A, 9 
in B). However, all the Fe n lines are weak and the majority of 
them are in the UV range, were the S/N ratio of our spectra drops 
dramatically. 

The ((3D) - IDmarcs) corrections for Fen are in the oppo- 
site direction of to those for Fe 1, making the ionisation imbal- 
ance worse. The different signs of the ((3D) - IDmarcs) correc- 
tions for neutral and ionised species are due to the different ioni- 
sation structure of the (3D) and ID models. We note that for the 
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Fig. 7. Synthetic spectral fits to the Li line in the co-added U VES 
spectrum of star A (top) and B (bottom). The observed spectra 
have been corrected for veiling (see Table [5J, so the lines appear 
with their intrinsic strength in each star. 



metal-poor subgiant HD 140283. IShchukina et aTl (120051) have 
performed NLTE computations for Fe i and Fe n, using a single 
snapshot of a hydrodynamical simulation, and found 3D-NLTE 
corrections of +0.6 for Fe i and +0.4 for Fe n. To the extent that 
these computations can be considered representative of the stars 
in CS 22876-032, we expect that a full 3D-NLTE analysis might 
achieve a better ionisation balance for iron. 

5.2.4. Lithium 

The high quality of our spectra allowed us to measure the Li 
doublet in both components of CS 22876-032 for the first time. 
Fig.[7]shows our mean spectra of the Li region, where our well- 
resolved spectra (only) have been co-aligned on the lines of star 
A (top) and B (bottom), respectively. The superposed lines of 
the other star appear only slightly diffuse, because the velocity 
difference between the two stars is nearly the same in all our 
UVES spectra (see Table [TJ. The spectra in Fig. [7] have been 
corrected for veiling, as discussed in § 15.2.11 

The Li abundance of each star was computed from two 
sets of spectra taken on different nights. We found differences 
of 0.01 and 0.05 dex between star A and B and adopted 
the average value of the two measurements. Table [5] lists the 
observed equivalent widths and the average Li abundances, 
with an d without the NLTE corrections obtained from the ta- 
bles of Carls son et al. I (Il994l) . We further correct for the ef- 
fect of depletion as predicted by t he sta ndard isochrones of 
iDelivannis. Demarque. & Kawaler I ([1990). The correction is 
negligible for the hot primary, but larger for the secondary (see 
TableU). 
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It is certainly surprising that the two components appear 
to have a different lithium content. The Li abundance of the 
primary component seems to be consistent with that observed 
in othe r metal-poor stars, i.e., the Spite plateau dSpite & S pite 
1982a|bT: iBonifacio & Molaro I 1 19971 iRvan. Norris. & B eers 
19991; iMelendez & RarmrezH2004t ICharbonnel & Primas Il2005t 
Asplund et al. 1 12006; Bon ifacio et al. I I2007I) . modulo the uncer- 



tainties on the temperature scales adopted over the years by the 
various authors, but the secondary star definitely seems to ex- 
hibit a lower Li content. 

According to the Li d e pletion isochrones of 
Deli vannis. De marque, & Kawalerl 0.990), if star B were 
350 K cooler than we assume, the correction for A(Li) would 
be 0.6 dex. Such a change in T e $ would also imply a slightly 
higher logg from the isochrone, so the model dependencies 
given in Table |4] would imply that A(Li) should be reduced 
by 0.17, giving a "corrected" Li abundance of A(Li)=2.20, 
in essential agreement with A(Li) of the primary (we do not 
consider variations in the temperature of star A). By changing 
simultaneously the effective temperatures of the pair (within 
the range allowed by photometry), we would derive different 
Li abundances for each of the components. However, we 
would also change the veiling factors, and the final abundance 
difference between the two stars would be roughly the same. 

We have checked whether inaccurate veiling corrections, es- 
pecially for the fainter lines of star B, could be the cause of 
the different Li abundances. It turns out that this is impossible: 
bringing A(Li)a into consistency with the Spite plateau would 
require a doubling of the veiling correction at 670 nm; given 
that we find consistent abundances for Na at 590 nm and Mg I 
at 880 nm, the veiling correction cannot be off by a factor two at 
the intermediate wavelength. 

For the Li lines we also performed a full 3D-LTE syn- 
thesis using Linfor3D . The resulting corrections are listed in 
Table [5] Since the 3D computation was performed in LTE, the 
Li abundances must no t be taken as defini t ive, a s shown by 
ICavrel & Steffenl {2000) and I Asplund et al. I d2003l) . It is, how- 
ever, interesting to notice that while the 3D effect in star A is 
almost entirely due to the horizontal temperature fluctuations, 
for star B it is almost entirely due to the cooler average tempera- 
ture profile of the 3D model. A full 3D-NLTE synthesis of Li in 
CS 22876-032 is beyond the scope of this paper. Full 3D-NLTE 
synthesis of the Li profile in HD74000 have been addressed re- 
cently bv lCavrel et ail d2007l) . 



5.2.5. Oxygen 

The oxygen abundances have been derived from UV OH 
lines of the (0-0) vibrational band of the A 2 E - X 2 Yl elec- 
tronic system. The use of these lines for oxy gen measure- 
ments in metal-poor stars was pioneered by iBessell et al. I 
( 1 99 if)- We were able to measure four lines in the primary 
and nine lines in the secondary. Followin g the extens i ve sur - 
veys of OH lines in me t al-poor stars by llsraelian et all £1998); 
iBoesgaard et al. I d 19991) ; llsraelian et al. I d2001l) . and the contro- 
versial fin ding of strongly increasing fO/F el with decreasing 
metallicitv. lAsplund & Garcia Perez 1 Q2001 ) warned of the pos- 
sible role of 3D effects on the formation of these lines. For this 
reason we decided to compute ad hoc CO s BOLD hydrodynam- 
ical simulations for this system in order to correctly evaluate the 
3D effects. Our results are summarised in Table [6] Some of the 
analysed lines are shown in Figs.[8ll9land[T0l 

The 3D effects are cl early large and considerably di ffer- 
ent between the two stars. lAsplund & Garcia Perez I d2001l) at- 




Fig. 8. Synthetic spectral fits to OH lines in the co-added U VES 
spectrum of star A (top) and B (bottom) in CS 22876-032. The 
observed spectra have been corrected for veiling factors of 1.30 
(A) and 4.36 (B). 




Fig. 9. Same as Fig. |8]for another OH line. 




Fig. 10. Synthetic spectral fits to OH lines in the co-added 
UVES spectrum of star B. The observed spectra have been cor- 
rected for a veiling factor of 4.36. 
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Table 6. 3D abundance corrections for the OH lines. [O/H] and 
[O/Fe] are relative to loge(O) = 8.72. 



A (nm) 


[O/H]^ 


[0/Fe] 1D A 


3D-<3D> 


A3D-1D 


CS 22876-032A 


311.21 


-1.48 


2.18 


-0.633 


-1.546 


311.47 


-1.55 


2.11 


-0.536 


-1.328 


312.26 


-1.56 


2.10 


-0.590 


-1.420 


313.43 


-1.48 


2.18 


-0.797 


-1.664 


CS 22876-032B 


312.26 


-1.63 


1.94 


0.027 


-0.985 


312.77 


-1.85 


1.72 


0.014 


-0.868 


312.83 


-1.68 


1.89 


0.033 


-1.009 


313.03 


-1.94 


1.63 


-0.004 


-1.045 


313.32 


-1.69 


1.88 


0.003 


-0.856 


313.66 


-1.75 


1.82 


-0.015 


-0.849 


314.30 


-1.80 


1.77 


-0.011 


-0.815 


317.45 


-1.75 


1.82 


-0.017 


-0.916 


317.77 


-1.70 


1.87 


-0.009 


-0.906 



Table 7. ID gravity effects on the Fe and O abundances 



Species 


Alogg 


A[X/HL 


A[X/H] S 


OH 


0.3 


-0.12 


-0.08 


Fc 1 


0.3 


-0.01 


-0.02 


Fen 


0.3 


0.09 


0.11 



tributed the 3D corrections primarily to the different average 
temperature profile of the 3D models, in particular to their ex- 
tremely cool outer layers. However, we see that for the primary 
star the 3D correction is almost evenly shared between average 
temperature profile and horizontal temperature fluctuations^. 

Our cooler 3D model atmosphere exhibits much smaller tem- 
perature fluctuations around the mean than the hotter one. In 
the cool component the average temperature becomes so low 
that a substantial amount of H2 is formed in the higher photo- 
spheric layers. The associated increase of the specific heat makes 
it much harder for pressure fluctuations to introduce tempera- 
ture fluctuations. Due to the smaller temperature fluctuations the 
cooler 3D and (3D) models provide essentially the same abun- 
dances. The model of the primary star is hotter on average, H2 
molecules are much less abundant, and temperature fluctuations 
are much more pronounced. Consequently, the resulting abun- 
dances differ between 3D and (3D) models. 

Fig. [TT] illustrates the situation. The 3D model exhibits a 
stronger cooling with respect to the IDlhd model at 5900 K than 
is the case at 6500 K. Over a wide pressure range the structure 
is almost adiabatic and passes through a region of substantial H2 
molecule formation, indicated by the rather high values of the 
specific heat, suppressing temperature fluctuations in that region. 
The almost adiabatic structure of the cool model also indicates 
that the convective overshooting is very efficient compared to ra- 
diative heating, and has driven the thermal structure into almost 
adiabatic equilibrium. 

4 -0.64 dex due to temperature fluctuations and -0.85 dex due to 
average temperature profile 




log 10 (P/dyne cm' 2 ) 

Fig. 11. Entropy profiles of the (3D) C0 5 BOLD model at- 
mospheres (solid lines), compared to IDlhd models (dashed) 
for stars B (5900 K) and A (6500 K). The labelled contours 
indicate the specific heat at constant pressure in units of 
10 8 erg g -1 Kr'(see text). 

The large difference in the behaviour of our two models 
warns us that to measure reliable abundances from OH one needs 
a grid of 3D models which is fairly dense in temperature, to cap- 
ture, for any metallicity, the T e ff at which H2 formation sets in. 

Finally, we note that the O abundance derived from OH lines 
is rather sensitive to the surface gravity, as can be seen from 
Table [7] An increase of gravity of 0.3 dex introduces a decrease 
in O abundance of about 0.1 dex. This gravity dependence is 
larger than that of Fe 1, however it is advisable to use Fe 1 to de- 
rive [O/Fe], since the gravity dependence of Fen is in the oppo- 
site direction. Note that the values reported in Table|7]have been 
estimated without revising the veiling factor when changing the 
surface gravity of the models, contrary to the error estimates de- 
scribed in § 15.2.21 

5.2.6. Other elements 

The Ben resonance doublet at 313.0 nm is within the wave- 
length range covered by our spectra. The S/N ratio in that re- 
gion is ~ 25, yet none of the Be lines is detected, which 
implies an upper limit for the Be abundance log(Be/H) < 
-13.0. As this is an order of magnitude higher than the 
Be abundance expected for t hese stars from the t rend of Be 
abundance with metallicitv dGilmore et alj 119921: lRvanetal.1 
1 19921: iMolaro et al.lll997t iBoesgaard et al.lll999t iPrimasI 120021: 
Boesgaar d & Novickill2006l) . this result is not of deep signifi- 
cance. 

For the other elements, abundances were determined directly 
from the EWs, except for Sc and Co for which we used spectrum 
synthesis to take hyperfine splitting (HFS ) into account. Fo r Co 
we used the A and B factors measured bv lPickeringl d 19961). For 
the Sc 11 361.3 nm line we used the A and B factors measured by 
Gangrskv et al. 1(120061) . In both cases we used the code linestruc 
of IWahlgrenl (12005) to compute HFS components. For the Sen 
42 4.6 nm fine we us e d the HFS components given in Table 5 
of iMc William et all d 19951) . The detailed atomic data of HFS 
for Sc and Co are provided in Table IA.2I For Mn we used a 
single line, the strongest of the Mn 11 lines of Mult. 3, for which 
hyperfi ne splitting is negligible according to Castelli & Hubrig 
(2004). 
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Silicon was only measured from the Si I line at 390.5 nm. In 
metal-poor cool giants, the silicon abundance is derived from a 
line at 410.3 n m, since the line Si I 390.5 nm is severely blended 
with CH lines dCavrel et al. 1120041) . The Si I at 410.3 nm is very 
weak in metal-poor dwarfs, but CH lines are so weak that the 
line at 390.5 nm can be used. 

Aluminium was measured from the resonance lines All 
394.4 nm and All 396.1 nm. The All 394.4 nm is blended with 
CH lines, which are also extremely weak in metal-poor dwarfs 
and were not taken into account. In fact, we computed synthetic 
spectra of the CH lines at 390.5 nm and 394.4 nm for the stellar 
parameters of both dwarfs, and the CH lines are not visible even 
for [C/Fe] as large as +2 dex. However, the 394.4 nm line pro- 
vided a significantly larger abundance than that derived from the 
396.1 nm line in both stars. For this reason, in Table [3] we give 
the average Al abundance derived from both lines with a large 
dispersion of ~ 0.2 dex. 

Using (3D) models, we determine the ((3D) - IDmarcs) cor- 
rections listed in Table[3] We note that the ((3D) - IDmarcs) cor- 
rections also take into account different veiling factors estimated 
using (3D) models and ID OSMARCS models. In general, the 
veiling factors estimated from (3D) models are higher for star 
B and lower for star A than those obtained using ID models. 
This effect is especially important for the secondary star, and be- 
comes more significant at shorter wavelengths. Taking into ac- 
count this effect, it is interesting to note that ((3D) - IDmarcs) 
corrections are in general negative for neutral species and pos- 
itive for ionised species, at least for the primary star, for which 
the veiling factors are not significantly different from those esti- 
mated using ID models, as already noted for iron. This reflects 
the different ionisation structure of ID and (3D) models. The 
((3D) - IDmarcs) corrections do not help to achieve ionisation 
equilibrium for Ca, Cr and Fe, the only elements for which abun- 
dances were determined from both neutral and ionised species. 

6. Discussion 

6.1. Lithium 

In Fig.[T2]we show the lithium abundances for the two compo- 
nents of CS 22876-0 32, together wit h our d ata from Paper VII, 
including the data of Asplun cTet al. I (|2006) rescaled in [Fe/H] 
and A(Li) to be homogeneous with our own data. Star A ap- 
pears to have a Li abundance at the same level as the major- 
ity of stars with metallicity below -2.5, if anything slightly 
higher. Star B appears to be far belowQ any of the other mea- 
sured stars. It should be noted that all the other stars in Fig. [12] 
have effective temperatures de termined from the wi ngs of Ha 
using the broadening theory of Barkler rTet al. I (|2000), while for 
CS 22876-032 they have been determined independently from 
colours and isochrones. However, the reasonable agreement be- 
tween our computed Ha profile and the observed profile shown 
in Fig.|6]suggests that the two temperature scales are fairly close. 

The fact that at the lowest observable metallicity, star A re- 
mains at the level of the plateau suggests that there is no down- 
turn or decrease in Li abundance at the lowest metallicities. This 
suggests that the slope of A(Li) with [F e/H] which is detectable 
in the sample of Asplund et al. (2006) (but not in that of Paper 
VII, alone) is not real, but rather an artifact due to the Ha tem- 
perature scale. It is possible that it is ultimately due to our inabil- 
ity to correctly model the atmospheres of extremely metal-poor 



5 -0.33 dex below the average Li abundance of the stars presented in 
paper VII 
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Fig. 12. 1D-NLTE Li abundances vs. [Fe/H] for the stars in 
CS 22876-032 (circles) and in other metal-poor dwarfs as 
reported by Asplund et al. (2006, triangles) and Bonifacio 
et al (2007, rhombs). The Li abundances by [Asplund et al. I 
(2006) were rec omputed using turbospectrum as reported by 
Boni facio et al. I d2007l) . No 3D-LTE or 3D-NLTE has been con- 
sidered. 
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Fig. 13. Same as Fig. |T2l but in the Li-r e ff plane. 

stars and the wings of Balmer line profiles. In Paper VII we ar- 
gued that the data could suggest either a vertical drop or an in- 
creased scatter in A(Li) at the lowest metallicities. Now, the drop 
seems to be ruled out by the A(Li) measured in star A. The mea- 
surement in star B, taken at face value, may support the idea that 
at metallicities below -2.5, the Spite plateau displays a sizeable 
scatter. 

In section |5~.2.4| we have pointed out how the difference of the 
Li abundance between the two stars could be resolved by assum- 
ing that the temperature of star B were 5550 K. Star B would thus 
be subject to considerable Li depletion, according to the standard 
Li dep letion isochrones of Delivanni s. Demarque. & Kawaler I 
( 1990). Such isochrones are available only for metallicities con- 
siderably higher than that of CS 22876-032. If, for any reason, 
either the lower metallicity of our system, or inclusion of other 
physical phenomena, the dependence of Li depletion on T e g is 
steeper than predicted by purely diffusive standard isochrones, 
then to reconcile the Li abundances of the two stars, the temper- 
ature of star B could be higher than 5550 K. Considering that 
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our estimated error on the effective temperature of star B is 150 
K {2(f), such cooler temperatures are not totally implausible. 

In our view the existence of a real scatter in Li abundances 
at the lowest metallicities remains to be established beyond any 
reasonable doubt. It is, nevertheless, worthwhile to discuss the 
possible implications of such a scatter, if real. 

We have already noted in Paper VII that by arbitrarily divid- 
ing the sample into two sub-samples, one with metallicity below 
-3.0 and the other above, the scatter of the "higher" metallicity 
sample was 0.05 dex, while the scatter of the "lower" metallicity 
sample was 0.11. If to the lower metallicity sample we add the 
two stars of CS 22876-032, th e scatter increases only slightly, to 
0. 12 dex. lRichard et al.l d2005l) have studied the effect of atomic 
diffusion in presence of turbulence and concluded that the ob- 
servations of the Spite plateau could be explained by starting 
from a primordial A(Li), compatible with the baryonic density 
derived from the WMAP experiment. The effect of atomic dif- 
fusion, countered by a suitably parametrised turbulence, can ex- 
plain the present level and low scatter of the Spite plateau. 

From their Li i sochrones in the pres ence of pure atomic dif- 
fusion (Fig. 5 of Rich ard et al. I l200"5l) it is obvious that one 
should expect Li abundances in the range 1.70<A(Li)<2.35. The 
sample constituted by the 8 stars from Paper VII with metal- 
licity below -3.0 and the two stars in CS 22876-032 spans the 
range 1.91 - 2.20 in A(Li). One could therefore suspect that in 
extremely metal-poor stars, turbulence is lower and atomic dif- 
fusion more efficient, causing the increased scatter in Li abun- 
dance. In Fig. Q~3] we see no clear tren d of A(Li) vs. r e ff. W hen 
comparing this figure with Fig. 5 of iRichard et al. I d2005l) we 
note that the "high" Li abundances of the hottest stars in the sam- 
ple (among which is the primary of CS 22876-032), preclude 
any clear resemblance between the two pictures. Therefore, cur- 
rent pure diffusion models seem unable to explain at the same 
time the behaviour of Li abundances with [Fe/H] and r e ff. It 
is possible, however, that they may do so, after an ad hoc 
parametrisat ion of turbu l ence with metallicity. 

Re cently iKorn et all (120061 l2007t) have found that the mod- 
els of IRichard et al. I d2005l) . with a suitable value for the tur- 
bulent diffusion coefficient, can explain the 0.12 dex difference 
in A(Li) they find between stars at the turn-off and on the sub- 
giant branch stars of the globular cluster NGC 6397. The same 
authors, however point out that assigning temperatures for the 
TO stars, hotter by 170 K (there fore close to the temperatures 
adopted by Bonifac io et alj200 2 for the TO stars of this cluster), 
such a difference would vanish. 

Though suggestive, the applicability of such turbulent diffu- 
sive models remains to be proven. The main cause of concern is 
the parametrisation of the turbulent diffusive coefficient, which 
is linked to a fi xed temperature, an d not to the bottom of the 
convective zone dRichard et al.ll2002h . 

In order to strengthen the observational constraints on such 
models and refine the estimates of the scatter and slope (or lack 
thereof) of the extreme metal-poor end of the Spite Plateau, fur- 
ther high-quality spectroscopy of EMP stars and additional accu- 
rate constraints on the effective temperatures of the whole sam- 
ple are highly desirable. 

6.2. Oxygen 

Fig. [l4]compares our 3D-LTE O abundances obtained from OH 
UV bands in CS 22876-032 with the high-quality measurements 
of Paper V for giant stars (mixed and unmixed), obtained from 
the [OI] 630nm line. We note that we have computed the 3D- 
LTE [O/Fe] ratios using the true 3D-LTE oxygen abundances 
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Fig. 14. 3D [O/Fe] ratios of the stars in CS 22876-032 (cir- 
cles) and in the metal-poor giants of Paper V. Triangles: "un- 
mixed", rhombs "mixed" stars; downward triangle with arrow: 
upper limit for CS 22172-002). 



and the (3D) Fe i abundances. We have chosen Fe i rather than 
Fe n abundances as the reference, because their sensitivity to the 
surface gravity is similar as that of the OH bands. The reader 
might wonder why we find different O abundances for stars A 
and B, but they are in fact consistent within the error bars, which 
mainly reflect the uncertainties in the effective temperature (see 
Tableg}. 

The abundances derived from [OI] are the raw values ob- 
tained from 1D-LTE model atmospheres. One worry is the pos- 
sible effects of granulation on the abu ndances derived from the 
[OI] line in giants. In a recent paper. ICollet et al.l d2007l) investi- 
gated the 3D effects in giant stars and found very small correc- 
tions for the [OI] lines for metallicities down to -2.0, but sizeable 
(almost 0.2 dex) downward corrections for [O/Fe]) at metallicity 
-3.0 (their Fig. 13 and Table 3). If we appl y the corrections in- 
terpolated and extrapolated from Table 3 of ICollet et alj {2007) 
to the giants of Paper V, the [O/Fe] decreases in all the giants 
and the good agreement between CS 22876-032 and the giants 
no longer holds. The mean [O/Fe] of the giants would be 0.51, 
while that of CS 22876-032 is 0.87. If instead we take the mea- 
surements of Paper V at face value, the mean [O/Fe] is 0.72. 

At present we have no full 3D models for giant stars, how- 
ever, from a few snapshots for models of r e jf = 4900, log g = 2.0 
and metallicity -2.0 and -3.0 we find little difference at the two 
metallicities, and very small 3D corrections. The 3D - IDlhd 
abundance correction amounts to -0.006 dex at metallicity -2.0 
and -0.037 dex at -3.0, which is negligible. For our CO s BOLD 
models of giants the difference between the mean 3D temper- 
ature structure and a corresponding IDlhd (<*mlt= L0) model 
is not large. In particular we do not find the substantial cool- 
ing of the highest p hotospheric layers at metallicity -3.0, which 
Coll et et all d20071) find (see their Fig. 1). Moreover, the mean 
3D temperature profile is slightly hotter than the IDlhd tem- 
perature profile in the [OI] 630 nm line-forming layers around 
log t ~ - 1 . The differences between our assessment of the im- 
pact of granulation on the abu ndances der i ved fr om the [OI] 630 
nm line in giants, and that o f ICollet et all d2007), will be further 
investigated in the future. This might be rooted in the different 
binning schemes adopted by the two codes for the opacity (4 
opacity bins for the Stein & Nordlund code and 6 opacity bins 
for CO s BOLD). In any case we believe that, at present, it is safer 
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not to apply any 3D correction to the O abundances derived for 
giants from the [OI] 630 nm line. 

In Paper V we found that the mean value of [O/Fe] in the 
range -3.6 < [Fe/H] < -2.5 was about 0.72 ± 0.17. The 
weighted average value in CS22876-032, [O/Fe] = 0.83 + 0.15, 
agrees with this determination. This mean value does not ex- 
clude a slight increase of the ratio [O/Fe] in the range -3.6 < 
[Fe/H] < -2.5 as seen in the theoreti cal predictions of th e 
chemical evolution models presented by Franc ois"et al.l (|2004). 
A more detailed discussion of the general behaviour of oxygen 
with metallicity will be treated in forthcoming investigations, 
hopefully after 3D corrections will be determined for all the 
stars, and for different lines of oxygen, Fe i and Fe n. 

6.3. The Odd-Z light elements 

The LTE [Na/Fe] abundance ratios of both stars appear con- 
siste nt with the Galacti c trend of this element in EMP gi- 
ants (ICavrel et al. Il2004|). subgiants an d turn-off stars of simi- 
lar metallicity dAndrievskv et al. Il2007l) . NLTE corrections are 
expected to be larger for the metal -poor giants than for dwarfs. 
We have estimated the NLTE correction^ for Na to be Anlte ~ 
-0.06 dex , according to th e NLT E corrections reported in 
Table 2 of lAndrievskv et aTl d2007l) . After applying these cor- 
rections, Na abundances of the CS 22876-032 dwarfs remain 
compatible with those of metal-poor giants and dwarfs, which 

exhibit almost a constant ratio [Na/Fe] 0.20 in the metallicity 

range -4.0 < [Fe/H] < -2.5. 

The LTE [Al/Fe] abundance ratios of the dwarfs in this sys- 
tem are ~ 0.5 dex larger than those in metal-poor giants with 
similar iron content. However, alu minium is also expected to ex - 
hibit significant NLTE corrections (Bau muller & Gehren 111 9971) . 
which might explain this difference, as was the case for Na. This 
element also exhibits an almost constant ratio [Al/Fe] ~ -0.10 
for giants in the metallicity range -4.0 < [Fe/H] < -2.5 when a 
fixed NLTE correction of +0.65 dex is considered. 

Within the errors, [Sc/Fe] is consistent with, although 0.2- 
0.3 dex lower than, the [Sc/Fe] ratio in EMP giants which shows 
an almost constant ratio [Sc/Fe] ~ 0. Note that the [X/Fe] ra- 
tios given in Table [3] were computed relative to Fe n for ionised 
species. As noted in § 15.2.31 the Fe n abundances are less reliable 
than those for Fe i because the Fe n lines in EMP dwarfs are very 
weak, especially for star B in CS 22876-032. Therefore, [X/Fe] 
for ionised species in star B should be regarded with caution. 

6.4. The a elements 

The [Mg/Fe] ratios in both stars of CS 22876-03 2 seems to be 
fairly consistent with those f ound in EMP giants (ICavrel et al. I 
2004), subgiants and dwarfs dCohen et al. 1 120041) . Bonifacio et 
al. 2007, in prep.), at the level of ~ 0.3 - 0.4 dex. 

[Ca/Fe] was derived from both Cai and Can lines, which 
yield similar [Ca/H] abundances, at least for the primary star. 
However, although the [Ca/Fe] ratios differ by ~ 0.2 dex, they 
seem to be slightly lower than those measured in metal-poor gi- 
ants and dwarfs, where a constant [Ca/Fe] ~ 0.4 is seen. 

[Si/Fe] also seems to be low, at [Si/Fe] 0.1 for stars A and 

— 1-0.1 for the secondary, compared to the constant [Si/Fe] ~ 0.4 
for metal-poor giants. However, other metal-poor dwarfs show 
similar abundances, which might be related to the different Si 
lines used in giants and dwarfs (Bonifacio et al. 2007, in prep.). 



6 A NLTE = log e(X) NLT E - log e(X) LTE 



Finally, [Ti/Fe] is constant at ~ 0.3 in EMP giants. While we 
find a similar result for star A, [Ti/Fe] in star B is completely dif- 
ferent due to its high Fe n abundance. We note that the standard 
deviation in [Ti/Fe] is relatively high in both stars (cr ~ 0.3), 
although the mean abundances are better determined when aver- 
aging the results from the 19 lines in star A and 12 in B. 

6.5. The iron-peak elements 

Chromium was derived from both Cr i and Cr n lines, and we find 
a difference of ~ 0.2 dex for the [Cr/Fe] ratios. [Cri/Fe i] in CS 
22876-032 appears to agree with that in other EMP dwarfs, but 
is slightly higher than seen in EMP giants. 

The ratio [Mn/Fe] is consistent with the other EMP giants 
and dwarfs, at [Mn/Fe]~ -0.5. [Ni/Fe] also agrees with the val- 
ues reported for EMP giants, which show a constant [Ni/Fe] ~ 0. 

Cobalt is found to be slightly enhanced in CS 22876-032 rel- 
ative to the gradually increasing trend of [Co/Fe] with decreas- 
ing [Fe/H] observed in giants, although marginally compatible 
within the errors. 

6.6. Could CS 22876-032 be a triple system ? 

One might ask whether a third star in CS 22876-032 might con- 
tribute significantly to the total light and the veiling of the lines 
of the two main components we have discussed so far. Such a 
star would need to have a mass above ~0.5 M in order to have 
any significant effect on the observed spectrum. The presence of 
such a third star can be ruled out by two independent pieces of 
evidence. 

First, we have individual spectra of CS 22876-032 with S/N 
ratios ~ 100 in the region of the Mg i b triplet. These are among 
the strongest stellar lines seen in these EMP stars and would be 
at least as strong in the third star. The Cayrel formula predicts 
that any line with an EW above ~ 0.15 pm would be detected 
at the 3<x level. Assuming teff = 5500K and logg = 4.7 (from 
the isochrone) plus [Fe/H]= -3.6 and [Mg/Fe] = +0.2 for the hy- 
pothetical third star, we need to dilute the strongest line of the 
Mgi b triplet at 518.36 nm which would show an intrinsic EW 
of ~ 15.0 pm. In Appendix lAl we define the veiling factors for a 
triple system and conclude that the non-detection of the 518.36 
nm line requires > 100. We can therefore conclude that any 
third star contributes negligibly to the total light of CS 22876- 
032. 

A second line of evidence is available from the radial ve- 
locities, which are accurate to ~1 km s _1 and cover a period of 
16 years. Any third star of mass comparable to A and B should 
leave significant trends in the velocity residuals from the orbital 
solution for periods shorter than several decades. We have there- 
fore attempted sinusoidal fits to the velocity residuals for both 
stars and find periods of the order of 1300 days in both cases. 
First, the standard deviations around these fits are 0.83 and 0.51 
km s _I , well below the purely observational errors, which shows 
that these results cannot be statistically significant. Second, a 
period ratio of ~ 3 between the outer and inner orbits is far too 
small for a triple system to be dynamically stable. A white dwarf 
in an orbit of much longer period is a possibility, but would not 
be detectable in our spectra. 

In summary, we conclude that the abundance results reported 
here cannot be significantly affected by light from a third star in 
the system - certainly not the discrepant Li abundances of the 
two stars. 
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7. Summary 

Our high-resolution VLT/UVES observations of the double- 
lined spectroscopic binary CS 22876-032 confirm that it har- 
bours the most metal-poor dwarfs known so far. Our improved 
orbital elements, together with published photometry and theo- 
retical isochrones, enable us to determine stellar parameters of 
T e ft; A = 6500 + 100 K and log g A = 4.4 + 0.1 for the primary 
(star A) and T eS ,B = 5900 + 150 K and log g B = 4.6 + 0.1 for 
the secondary (star B). 

Using ID OSMARCS models and the turbospectrum code, 
we determine abundances of Li, O, Na, Mg, Al, Si, Ca, Sc, 
Ti, Cr, Mn, Fe, Co, and Ni, correcting the observed spectra for 
the veiling from the continuum flux of the other star. We find 
[Fe/H]= -3.66 + 0.16 and [Fe/H]= -3.57 + 0.25 for star A and 
B, respectively. Using C0 5 BOLD model atmospheres to esti- 
mate 3D abundance corrections, we compute full 3D spectrum 
synthesis using the Linfor3D code for Li and O to estimate the 
3D - IDlhd corrections, while we use a horizontal and temporal 
average of the 3D model to compute (3D) abundances with tur- 
bospectrum for the rest of elements. In general, we find ((3D) 
- IDmarcs) corrections to be S> 0.1 dex, negative for neutral 
species and positive for ionised species; for Fe in particular, we 
find corrections to be ~ -0.12 (A)and 0.07 (B). 

The fa/ Fe] ratios are cons istent with our earlier results for 
EMP giants dCavrel et al. 1120041 First Stars V), although Ca and 
Si are slightly low ([X/Fe]$ 0), but actually consistent with 
our results for other EMP dwarfs (Bonifacio et al. 2007, in 
prep.). [Na/Fe] appears consistent with both EMP giants and 
dwarfs when NLTE corrections are considered. The LTE value 
of [Al/Fe] is not consistent with those in EMP giants, but Al 
is severely affected by NLTE effects, which may solve this dis- 
crepancy. The iron-peak elements follow the established trends 
in EMP giants and dwarfs. 

Our high-quality spectra allowed us to measure the Li dou- 
blet in both stars of CS 22876-032 for the first time. We find 
NLTE Li abundances of 2.18 + 0.05 and 1.77 + 0.09 for stars 
A and B, respectively. While the Li abundance of star A corre- 
sponds to the level of the Spite plateau, the secondary star has 
a significantly lower abundance. This discrepancy may be re- 
solved by assuming that the secondary star has been subject to 
significant Li depletion, which, according to standard Li deple- 
tion isochrones, would have been the case if the star were 350K 
cooler than assumed by our analysis. Full 3D corrections for Li 

are estimated to be 0.3 (A) and 0.2 (B); however, these 

computations were performed in LTE, and 3D NLTE corrections 
are needed to confirm the sign and value of these corrections. 

The near-UV part of our VLT/UVES spectra enabled us to 
measure oxygen abundances from the OH bands. We find ID 
[O/Fe] values of 2.14 + 0.15 (A) and 1.81 + 0.33 (B) and com- 
pute full 3D corrections for the OH lines, which turn out to be 
-1.5 (A) and -1.0 dex (B). Using these corrections and the (3D) 
Fei abundances, we determine 3D [O/Fe] ratios of 0.77 + 0.15 
(A) and 0.97 + 0.33 (B). These 3D [O/Fe] ratios are consistent 
with those derived from the [OI] line in EMP giants of similar 
metallicity, where 3D corrections should not be significant. 
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A wavelength 

the flux spectrum s A of 
the i,k component 
the continuum flux spectrum c\ of 
the /,/, component 
the line depression = c,- - s,- 
the intrinsic equivalent width of a 
spectral line of the /,/, component 
the observed equivalent width of a 
spectral line of the /,/, component 
the observed equivalent width of the 
three components. 
In any orbital phase we have: 



Si (i =1,2,3) 

cS= 1,2,3) 

di(i =1,2,4) 
EWi(i= 1,2,3) 

EW obs (i =1,2,3) 



EW obs 



EW 



pA2 
J/11 



d\ + d2 + di 
Cl + c 2 + c 3 



dA 



Where the interval [Al,<i2] includes the desired spectral 
feauture, which we assume to be isolated. Considering that 
the continuum fluxes can be assumed to be constant over the 
integration interval this can be re-written as: 



ew: 



obs 
123 



JAl Cl C]+c-2+C 3 XI 1 c 2 



.12 . 



JAl c 3 



We now define the veiling factors as f, = (ci + c 2 + c 3 )/q 
and consider a phase in which the radial velocities are such that 
the line of of each component is isolated and not blended with 
the others, then one has: 



EW 



obs 
' 123 

and noting that : 



EW, 



obs 



■ EWf s + EWf s 



r>A2 
J/11 



*dA 

Ci 



EW, 



one has that : 



EWf s + EW obs + EW° bs 



EW l /f 1 +EW 2 /f 2 + EW 3 /f 3 



finally considering that in this phase the lines are not over- 
lapping this equality implies: 



EWi = fiEW\ 



obs 



Appendix A: Veiling factors for a triple system. 

In this appendix we derive the expression for the veiling factors 
of a triple system. In fact they are a trivial extension of those for 
a double system, however since they are not readily found in any 
paper or book we know of, we provide them here for the reader's 
convienience. 

We shall use the following notation: 
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Table A.l. Line data, equivalent widths, veiling factors and ID abundances of CS 22876-032A,B 



Specie 


A 


X 




EW A ,obs 


/A, ID 


A(X) A , 1D 


EW B ,obs 


/b,id 


A(X) b ,id 




(nm) 


(eV) 




(pm) 






(dex) 


(pm) 




(dex) 


Nai 


588.9951 


0.000 


0.112 


1.40 


1 


.35 


2.63 


— 


— 


— 


Nai 


589.5924 


0.000 


-0.191 


0.79 


1 


.35 


2.64 


0.49 


3.87 


2.54 


Mgi 


333.6674 


2.720 


-1.230 


1.07 


1 


.32 


4.57 


0.90 


4.16 


4.68 


Mgi 


382.9355 


2.710 


-0.207 


5.01 


1 


.31 


4.41 


2.06 


4.24 


4.29 


Mgi 


383.2304 


2.710 


0.146 


6.03 


1 


.31 


4.29 


2.87 


4.25 


4.35 


Mgi 


383.8290 


2.720 


0.415 


6.89 


1 


.30 


4.21 


2.77 


4.31 


4.06 


Mgi 


416.7271 


4.340 


-1.000 


0.23 


1 


.29 


4.84 


0.11 


4.41 


4.86 


Mgi 


517.2684 


2.710 


-0.380 


4.88 


1 


.33 


4.47 


2.42 


4.04 


4.41 


Mgi 


518.3604 


2.720 


-0.158 


5.73 


1 


.33 


4.43 


3.24 


4.03 


4.52 


Mgi 


552.8405 


4.340 


-0.341 


0.60 


1 


.34 


4.58 


0.16 


3.94 


4.25 


Mgi 


880.6756 


4.340 


-0.137 


0.81 


1 


.39 


4.41 


0.62 


3.57 


4.50 


Ah 


394.4006 


0.000 


-0.623 


1.51 


1 


.28 


2.74 


1.12 


4.54 


2.88 


Ah 


396.1520 


0.010 


-0.323 


1.41 


1 


.31 


2.42 


1.16 


4.24 


2.56 


Si i 


390.5523 


1.910 


-1.090 


1.55 


1 


.29 


3.80 


1.32 


4.42 


4.07 


Cai 


422.6728 


0.000 


0.240 


3.60 


1 


.30 


2.71 


1.70 


4.39 


2.68 


Can 


317.9331 


3.150 


0.512 


4.08 


1 


.30 


2.62 


1.36 


4.28 


2.57 


Can 


370.6024 


3.120 


-0.480 


1.41 


1 


.34 


2.58 


— 


— 


— 


Can 


373.6902 


3.150 


-0.173 


3.22 


1 


.33 


2.84 


0.68 


4.04 


2.53 


Sen 


361.3829 


0.020 


0.416 


0.42 


1 


.33 


-0.69 


0.42 


4.01 


-0.40 


Sen 


424.6822 


0.310 


0.240 


0.37 


1 


.30 


-0.44 


0.16 


4.38 


-0.41 


Tin 


316.8532 


0.150 


-0.310 


2.35 


1 


.30 


1.66 


0.70 


4.29 


1.34 


Tin 


323.4520 


0.050 


0.426 


4.55 


1 


.31 


1.54 


1.64 


4.25 


1.34 


Tin 


323.6578 


0.030 


0.234 


4.11 


1 


.31 


1.56 


1.72 


4.25 


1.59 


Tin 


324.1994 


0.000 


-0.045 


3.20 


1 


.31 


1.51 


1.75 


4.24 


1.86 


Tin 


325.1918 


0.010 


-0.579 


1.54 


1 


.31 


1.52 


0.65 


4.23 


1.39 


Tin 


332.2941 


0.150 


-0.093 


2.38 


1 


.32 


1.44 


1.23 


4.17 


1.53 


Tin 


332.9453 


0.140 


-0.274 


2.53 


1 


.32 


1.65 


1.33 


4.17 


1.78 


Tin 


338.0279 


0.050 


-0.570 


3.38 


1 


.32 


2.12 


1.33 


4.14 


1.98 


Tin 


338.7846 


0.030 


-0.432 


2.29 


1 


.32 


1.63 


1.01 


4.14 


1.54 


Tin 


344.4314 


0.150 


-0.810 


1.02 


1 


.32 


1.63 


— 


— 


— 


Tin 


345.6388 


2.060 


-0.230 


0.41 


1 


.32 


2.32 








Tin 


347.7187 


0.120 


-0.967 


1.17 


1 


.32 


1.83 


0.58 


4.08 


1.76 


Tin 


348.9741 


0.140 


-1.920 


0.40 


1 


.33 


2.27 








Tin 


350.0340 


0.120 


-2.020 


0.37 


1 


.33 


2.32 


— 


— 


— 


Tin 


375.9296 


0.610 


0.270 


2.39 


1 


.32 


1.33 


0.94 


4.16 


1.26 


Tin 


376.1323 


0.570 


0.170 


2.36 


1 


.32 


1.39 


1.07 


4.16 


1.43 


Tin 


391.3468 


1.120 


-0.410 


0.39 


1 


.28 


1.48 








Tin 


402.8343 


1.890 


-0.990 


0.15 


1 


.29 


2.30 








Tin 


518.8680 


1.580 


-1.050 


0.14 


1 


.33 


2.01 


- 


- 


- 


Cri 


357.8684 


0.000 


0.409 


0.85 


1 


.ii 


1.79 


0.61 


4.03 


1.61 


Cri 


425.4332 


0.000 


-0.110 


0.45 


1 


.30 


1.83 


0.28 


4.37 


1.63 


Cri 


427.4796 


0.000 


-0.230 


0.32 


1 


.30 


1.79 


0.16 


4.36 


1.48 


Cri 


428.9716 


0.000 


-0.360 


0.46 


1 


.30 


2.09 


0.58 


4.35 


2.27 


Cri 


520.6038 


0.940 


0.020 


0.13 


1 


.33 


1.92 








Cm 


313.2053 


2.480 


0.451 


2.85 


1 


.30 


2.29 


1.10 


4.32 


2.37 


Mn ii 


344.1988 


1.780 


-0.270 


0.54 


1 


.32 


1.38 


0.32 


4.11 


1.53 
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Table A.l. Continued. 



opecie 


A 

(nm) 


X 




PW 

(pm) 


/a,id 


A (. A M,1D 
^uex ^ 


PW 

(pm) 


f 

JB,1D 


^uex ^ 


C Q t 

re i 


1AH C/1 CA 

j4/.j4jU 


A AQA 


1 AC/1 
-1.UJ4 


1 o/i 

J.Z4 


1 

1 


. JZ 


1 QQ 

j.yy 


1 OQ 

i.oy 


/I AQ 

4.uy 


/I AO 

4.U/ 


Fe i 


iah f^nt\^i 
J4/.0 /UZ 


A 1 OA 

U.1ZU 


1 CA"7 

-I.jU/ 


1 QA 

i.yo 


1 

1 


^0 
. JZ 


A AC 

4.Uj 


1 A A 
1.4U 


/I AO 
4.U0 


A AO 

4.U/ 


Fe i 


1 A o ACT/1 

j4y.Uj /4 


A ACA 

U.UjU 


1 1 AC 

-I.IUj 


O Q 1 

z.y 1 


1 
1 


.J J 


1 QA 

j.yu 


1 Q1 

i.yj 


A AO" 

4.U/ 


/i 1 1 
4.1 1 


P<3 T 

re i 


tci; CIOQ 

jjo.jj /y 


A QAA 

u.you 


A 1 11 
-U. 1 JJ 


O QA 

z.yo 


1 

1 


1,1, 
.J J 


1 7< 

J. /j 


1 QA 

i.yu 


A A/1 
4.U4 


1 QO 

j.yz 


T ; 

re i 


jjo. 1 iyj 


A OAA 

U.oOU 


A A 

U.4U0 


A ~\H 
4.1 / 


1 

1 


11 
.J J 


ICC 
J.JJ 


001 

Z.ol 


/I Al 
4.UJ 


1 OQ 

j.oy 


P^» T 

re i 


1AA OOCQ 

jou.oojy 


1 A1 A 
1.U1U 


A 1 AA 
-U. 1UU 


j.jy 


1 

1 


1,1, 
.J J 


1 Q 1 

j.y 1 


019 
Z.l / 


A A1 
4.U1 


A 1 1 
4. 1 J 


P^ T 

re i 


1A1 OOAO 
JOI.O /OS 


A QQA 

u.yyu 


a nm 
-U.UUj 


cc 
J.JJ 


1 
1 


11 
.J J 


1 oc 

J.OJ 


oc 
Z.Zj 


/I AA 
4.UU 


A AA 
4.U0 


T ; 

re i 


1A/1 -7O/II 

304. /o4j 


A QOA 

u.yzu 


A 1 Q/1 

-u. iy4 


1 OQ 

j.zy 


1 
1 


.J J 


1 -70 

5.11 


O O/I 

Z.Z4 


A AA 

4.UU 


A AQ 

4.uy 


P<3 T 

re i 


1 Q 1 CO/I A 

J01.J04U 


1 A QA 


U.Zj / 


7Q 

j. /y 


1 

1 


1,(1 


1 00 

J.OO 


1 A 1 
1.01 


A 1/1 
4. J4 


1 09 
5.11 


T ; 

re i 


oon A/1 OC 

J0Z.U4ZJ 


A OAA 

U.oOU 


A 1 1 Q 

u. 1 iy 


c 1 c 
J.l J 


1 

1 


. JU 


1 OQ 

j.oy 


O AA 

z.uu 


/I 1 A 

4.jU 


1 Al 

j.Oj 


T ; 

re i 


TO COO 1 

JOZ.JOOl 


A QOA 

u.yzu 


A A^T 

-U.Uj / 


/I 1 A 

4.1U 


1 
1 


1 1 
.j 1 


1 -7A 
J. 10 


O AA 
Z.UU 


A OA 

4.Z0 


1 OO 

j.oZ 


Fe i 


100 700'! 

JOZ. /OZj 


1 CAA 
l.JOU 


A AAO 
U.UOz 


O CO" 

Z.J / 


1 

1 


1,^ 
. j 1 


1 7Q 

j. /y 


1 OQ 

1 .zy 


/I oc 
4.ZJ 


191 
J. / 1 


C Q t 

re i 


10/1 A/110 
J04.U4J0 


A QQA 

u.yyu 


A CAA 
-U.jUO 


O OA 

z.zo 


1 

1 


. JU 


1 9A 
J. 10 


111 
1.11 


A 1/1 
4. J4 


1 AA 
J.OU 


T ; 

re i 


10/1 1 A/1 

Jo4. lU4o 


1 A 1 A 


A A/1 C 

-U.U4j 


O AC 

Z.Uj 


1 
1 


. JU 


1 -70 
J. /o 


A OO 
U.OO 


/lie 
4.JJ 


1 C/1 

J.j4 


re i 


10/1 qqao 

jo4.yyo / 


1 A1 A 
1.U1U 


A 57 1 

-U.o / 1 


1 OA 

1 .zo 


1 

1 


.zy 


1 9Q 

j. /y 


1 AQ 

i.uy 


A A A 
4.40 


1 QQ 

j.yy 


P^» T 

re i 


IOC AIOO 
JOJ.Oj /Z 


A ACA 
U.UjU 


1 OOA 
-l.ZOD 


1 O'Z 
J.Zj 


1 

1 


.Zo 


1 QA 

j.yo 


1 9A 
1./0 


/I CO 
4.JZ 


A 1 A 
4. 10 


Fe i 


Tor 001 1 

joj.yy i 1 


A AAA 
U.UUU 


A T 1 A 

-U. / 1U 


C A A 

J. 44 


1 

1 


98 
.ZO 


A Al 

4.Uj 


/1 
Z.4Z 


/I C/1 

4.J4 


A AO 
4.U0 


P^» T 

re i 


10A CCOI 

joO.jjZj 


1 Al A 
1.U1U 


A QOO 

-u.ysz 


1 A^ 
l.UJ 


1 

1 


98 
.ZO 


1 7Q 

j. /y 








P^ T 

re i 


07 OCA1 
jo /.ZjUl 


A QQA 

u.yyu 


A QOO 

-u.yzs 


1 CO 
l.JO 


1 

1 


9Q 
-Zv 


1 QC 

j.yj 


A QQ 

u.yy 


A /1Q 
4.4y 


1 Q/1 

j.y4 


T ; 

re i 


OA1 

Jo /.oUlo 


A OAA 

u.you 


A Q 1 A 

-u.y 14 


1 CQ 

i.jy 


1 
1 


. JU 


1 Q 1 

j.y 1 


1 OQ 

i.zy 


/I 1A 

4.J0 


AM 
4. 1 J 


r Q T 

re i 


100 AOOO 
JOO.OZOZ 


A ACA 
U.UjU 


1 ATA 

-1.U/0 


O A/1 

Z.04 


1 

1 


^9 

. JZ 


1 AA 
J.OU 


010 

Z.lo 


A 1 9 
4.1 / 


A 1 9 
4.1 / 


Fe i 


100 OA/1 Q 

JOO. /U4o 


A QOA 

u.yzu 


1 1/1/1 
-1.1 44 


1 AO 

1 .uz 


1 

1 


^9 
. JZ 


1 57 
J.O / 


11c 
l.lj 


/I 1 A 
4. 10 


/I 1 c 
4. 1 J 


T ; 

re i 


Ton CACA 

joy.jOjO 


A 1 1 A 

U.l 1U 


1 ATA 

-1.0/U 


1 CQ 

i.jy 


1 

1 


1 1 

. J 1 


1 QO 

j.yz 


1/11 
1.4j 


/I 1 Q 

4.iy 


/I 1 c 
4. 1 j 


P^» T 

re i 


1QQ Q7A7 

joy.y /u/ 


A AQA 

u.uyu 


1 C2 1 
-l.JJl 


O 1 A 

z. 10 


1 

1 


1,^ 
. j 1 


1 Q/1 

j.y4 


1 AA 
1.00 


/I OO 
4.Z0 


/I 91 
4.ZJ 


Fe i 


iqa OQ/1A 

jyu.zy4o 


1 CAA 
l.JOU 


A A 

-U.400 


1 1 A 
1.1U 


1 
1 


in 

. JU 


1 7Q 

j. /y 


A Q 1 

u.y 1 


A 1C 
4.JJ 


1 Q1 

j.yj 


Fe i 


IQA A/1 OA 

jyu.04oU 


A 1 1 A 
U.l 1U 


O/i Q 
-Z.Z4J 


A -7-7 
U.l / 


1 
1 


90 
.zy 


A AQ 

4.uy 


A Al 

U.Oj 


AAA 

4.44 


A AC 

4.UJ 


P^» T 

re i 


jyz.uzjo 


A 1 OA 
U.l ZU 


1 T/1 A 

-1. /40 


1 'JQ 

i.jy 


1 

1 


98 

.zo 


1 Q 1 

j.y 1 


1 OO 

1 .zo 


/ICO 
4. Jo 


A 9A 
4.ZU 


re i 


100 7oon 

jyz. /yzu 


A 1 1 A 
U.l 1U 


1 coo 
-l.jZZ 


O OA 

Z.ZO 


1 


.28 


1 fi7 

j.y / 


1 A A 

1.44 


/ICO 

4. Jo 


/I 1 
4. 1Z 


T ; 

re i 


/I A/1 CO 1 O 

4U4.J01Z 


1 A QA 

i.4yu 


A OOA 

U.ZsU 


^ QA 

j.yo 


1 


.29 


IOC 
J.OJ 


1 10 
l.JO 


a An 
4.4/ 


1 CO 

j.jZ 


Po T 

re i 


/1AA ICQ/I 
4U0. Jjy4 


1 CAA 

l.JOU 


A AAO 

u.uoz 


00 

Z.00 


1 

1 




1 0/1 
J.04 


1 A Q 

1 .4y 


/i /i 9 
4.4 / 


1 QA 

j.yu 


Fei 


407.1738 


1.610 


-0.022 


2.54 


1 

1 




3.87 


1.26 


4.44 


3.83 


Fei 


413.2058 


1.610 


-0.675 


0.98 


1 

1 


.zy 


3.96 


0.67 


4.42 


3.96 


P<3 T 

re i 


A 1 A 'JG/^Q 
4140o0o 


1 <AA 

l.JOU 


A C 1 1 
-U.Jl 1 


1 o/i 

1 .Z4 


1 

1 




1 57 
J.O / 


1 A9 
l.U/ 


A Al 
4.4J 


A AQ 

4.uy 


P^ T 

re i 


/ioa onoo 
4ZU.ZUZy 


1 A QA 


A TAO 

-U. /Uo 


A Ol 

u.yj 


1 

1 


9Q 
>Ly 


1 OC 
J.OJ 


A 7H 
U. /U 


A A A 
4.4U 


1 OQ 

j.oy 


Fe i 


/I 0< f\HQH 

4Zj.U/o / 


1 <aa 


A 11 A 

-U. / 14 


A CO 
U.JO 


1 
1 


JU 


1 AQ 

j.oy 


A /I Q 

u.4y 


/I 10 
4.J / 


1 O/I 

J. /4 


p^ t 
re i 


4Z0.U4 /4 


/iaa 
Z.4UU 


A 1 AQ 

u. iuy 


A TO 
U. /Z 


1 

1 




191 
J. / 1 


A OO 
U.Zo 


A 19 
4.J / 


i/ii 
J.4j 


C Q t 

re i 


a on 1 7 A 1 
4Z / . 1 /01 


1 A QA 


A 1 A/1 
-U. 104 


O /I ^ 
Z.4J 


1 

1 


JU 


1 OA 
J.OO 


1 01 
l.Zj 


A 1A 
4.J0 


1 OO 
J. 10 


T ; 

re i 


43Z.J /oZ 


l.OlU 


A AAA 

u.uuo 


O 1 Q 

z.iy 


1 
1 




1 71 

j. 15 


1 1Q 

i.jy 


/I 1 Q 

4.iy 


101 
J.Ol 


Fe n 


Jlo.O /3o 


1 o^a 
1. /UU 


1 7in 
-1 . / 1U 


1 AQ 

i.oy 


1 

1 


JU 


A 1 1 
4.1 J 


A QA 

u.yu 


A OO 
4.ZO 


A 19 
4. J / 


Fe ii 


Tin oqaq 

iiy.zyuy 


1.0 /U 


1 QCA 

-i.yju 


1 17 

Y.5I 


1 
1 


j 1 


A OA 

4.ZU 


A OQ 

u.zy 


/I 17 

4.Z/ 


1 OO 
J.O / 


Fe ii 


jly.j lyy 


1 Tin 
1. /ZU 


1 TOA 

-1. /ZU 


1 in 
I.jU 


1 


.31 


A A1 
4.U1 


A 7H 
U. /U 


/i 09 
4.Z/ 


A O 1 
4.Z1 


Fe ii 


J1V.OU/U 


1.0 /U 


1 AAA 
-1.00U 


1 A 1 
1.01 


1 


.31 


A Al 
4.UJ 


A A A 
U.44 


A 99 
4.Z / 


1 OA 
J.OU 


Fe ii 


101 (\A A A 

jZl. U444 


1 toa 
1. /ZU 


1 TQA 

-i. /yu 


O AC 

Z.Uj 


1 

1 


^ 1 


A 1 A 

4. jO 


A QC 

u.yj 


/I OA 

4.Z0 


A CO 

4.jZ 


Fen 


321.3309 


1.700 


-1.230 


1.62 


1 

1 


,D 1 


3.62 








Fen 


322.7742 


1.670 


-1.130 


3.29 


1 


.31 


4.05 


1.94 


4.26 


4.81 


Fe ii 


jZj.joo / 


a ooA 


-Z.jUU 


1 /IT 
1.4/ 


1 


.31 


A 1 Q 

4. iy 


U. /O 


A OO 

4.ZZ 


/I 1 1 

4.31 


Fe ii 


-5Z /. /34o 


A QQA 


-Z.4 IK) 


1 


1 


.31 


A 1 1 

4. 1 1 


A 70 
U. /Z 


A OA 

4.ZU 


/I O'* 

4.Z3 


Fen 


423.3172 


2.580 


-1.900 


0.51 


1 


.30 


4.26 








Fen 


492.3927 


2.890 


-1.320 


0.36 


1 


.32 


3.77 








Fen 


501.8440 


2.890 


-1.220 


0.42 


1 


.33 


3.73 


0.27 


4.08 


4.05 


Fen 


516.9033 


2.890 


-0.870 


0.65 


1 


.33 


3.58 








Co i 


340.5114 


0.430 


0.250 


1.17 


1 


.32 


1.93 


2.44 


4.12 


2.10 


Co i 


341.2333 


0.510 


0.030 


1.09 


1 


.32 


2.05 


0.86 


4.12 


2.02 


Co i 


345.3508 


0.430 


0.380 


1.54 


1 


.32 


1.99 


1.41 


4.11 


1.95 
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Table A.l. Continued. 



Specie 


A 


X 


log gf 


EW Aj0bs 


/a,\d 


A(X) A , 1D 


EW B ,obs 


/b,id 


A(X) B>1D 




(nm) 






(pm) 






^uex ) 


(pm) 




^uex ) 


Co i 


7. A O <£0 1 


a ^in 


-U.Z /U 


A 


1 


33 


Z.l 1 








CO I 


7CA 7770 

3jU.ZZ/o 


A A 7A 
U.43U 


a mn 


1 A A 
1 .4U 


1 


33 


7 11 
Z. 1 1 








Co i 




a non 

u.yzu 


A 77A 


A 7/1 

U.Z4 


1 


29 


7 AA 

z.uu 


A 1 O 

u.iy 


A A Q 

4.4o 


1 O/l 

l.y4 


Co i 


/i 1 1 17.11 

41Z. 131 1 


A o.7A 

u.yzu 


A 77A 
-U. 3ZU 


A 1 

U.lo 


1 


31 


1 OA 

i.yu 








INl I 


old OOOI 


a A7A 


A </1 A 

-U.34U 


7 


1 


32 


7 OA 

z.yu 


1 C7 


4.13 


7 07 

z.yz 


Mi t 

INI I 


7/1 7 ircyt 

343.33J4 


a 


-U.ooo 


1 CO 

i.sy 


1 


32 


7 77 
Z. / / 


1 A C 
1 .48 


A 1 A 
4. 1U 


7 Q7 

z.y3 


tm; r 
INI I 


7/1 C TQOC 

343.Z00J 


a 1 1 A 
U.l 1U 


A Q 1 A 

-u.y iu 


11/1 
1 . 14 


1 


32 


7 7Q 
Z. 10 


1.13 


/I 1 1 
4. 1 1 


7 Q7 

z.y3 


INl I 


34j.o4jO 


U.Z1U 


A Til 


1 7< 
Z.Zj 


1 


32 


7 /^7 
Z.OZ 


1 77 

1.3Z 


/I 1 A 

4.1U 


Z.M 


INl I 


34o.l64y 


A A7A 

0.030 


a 7 An 
-0.34/ 


7 c c 
Z.55 


1 


32 


7 £7 
Z.O / 


1 CI 

1.51 


/I 1 A 

4.10 


7 an 
Z.00 


Nil 


349.2954 


0.110 


-0.250 


2.39 


1 


33 


2.59 


1.76 


4.07 


2.84 


Nil 


351.5049 


0.110 


-0.211 


2.74 


1 


33 


2.67 


1.74 


4.06 


2.77 


Nil 


352.4535 


0.030 


0.008 


3.90 


1 


33 


2.77 


2.12 


4.05 


2.80 


Nil 


361.0461 


0.110 


-1.149 


0.95 


1 


33 


2.92 


1.07 


4.01 


3.05 


Nil 


361.9386 


0.420 


0.035 


2.02 


1 


33 


2.46 


1.72 


4.00 


2.75 


Nil 


380.7138 


0.420 


-1.180 


1.11 


1 


30 


3.17 


0.30 


4.38 


2.58 


Nil 


385.8292 


0.420 


-0.970 


0.73 


1 


28 


2.74 


0.68 


4.54 


2.85 
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Table A.2. Hyperfine structure of Sc and Co 



Specie A (nm) x (eV) logg/ Specie J(nm) ^ (eV) logg/ Specie i(nm) ;f (eV) logg/ 



Sc ii 


361.3815 


0.022 


-0.126 


Co i 


341.2328 


0.514 


-0.171 


Co i 


350.2241 


0.432 


-0.221 


Sen 


361.3814 


0.022 


-0.141 


Co i 


341.2325 


0.514 


-0.151 


Co i 


350.2263 


0.432 


-0.110 


Sen 


361.3819 


0.022 


-0.174 


Co i 


341.2333 


0.514 


-0.241 


Co i 


350.2255 


0.432 


-0.132 


Sen 


361.3817 


0.022 


-0.108 


Co i 


341.2330 


0.514 


-0.151 


Co i 


350.2245 


0.432 


-0.231 


Sen 


361.3815 


0.022 


-0.107 


Co i 


341.2325 


0.514 


-0.131 


Co i 


350.2272 


0.432 


-0.095 


Sen 


361.3823 


0.022 


-0.163 


Co i 


341.2338 


0.514 


-0.225 


Co i 


350.2262 


0.432 


-0.133 


Sc ii 


361.3821 


0.022 


-0.096 


Co i 


341.2333 


0.514 


-0.141 


Co i 


350.2251 


0.432 


-0.251 


Sc ii 


361.3817 


0.022 


-0.083 


Co i 


341.2326 


0.514 


-0.114 


Co i 


350.2283 


0.432 


-0.082 


Sc ii 


361.3828 


0.022 


-0.165 


Co i 


341.2343 


0.514 


-0.225 


Co i 


350.2271 


0.432 


-0.140 


Sc ii 


361.3825 


0.022 


-0.090 


Co i 


341.2336 


0.514 


-0.136 


Co i 


350.2257 


0.432 


-0.291 


Sen 


361.3821 


0.022 


-0.065 


Co i 


341.2328 


0.514 


-0.099 


Co i 


350.2295 


0.432 


-0.071 


Sen 


361.3835 


0.022 


-0.176 


Co i 


341.2349 


0.514 


-0.235 


Co i 


350.2281 


0.432 


-0.161 


Sen 


361.3831 


0.022 


-0.089 


Co i 


341.2341 


0.514 


-0.137 


Co i 


350.2308 


0.432 


-0.060 


Sen 


361.3826 


0.022 


-0.049 


Co i 


341.2330 


0.514 


-0.086 










Sc ii 


361.3843 


0.022 


-0.198 


Co i 


341.2356 


0.514 


-0.256 


Co i 


384.5473 


0.923 


-0.180 


Sen 


361.3838 


0.022 


-0.096 


Co i 


341.2346 


0.514 


-0.144 


Co i 


384.5475 


0.923 


-0.173 


Sen 


361.3832 


0.022 


-0.036 


Co i 


341.2334 


0.514 


-0.075 


Co i 


384.5470 


0.923 


-0.153 


Sen 


361.3853 


0.022 


-0.238 


Co i 


341.2364 


0.514 


-0.295 


Co i 


384.5479 


0.923 


-0.243 


Sen 


361.3847 


0.022 


-0.117 


Co i 


341.2352 


0.514 


-0.165 


Co i 


384.5474 


0.923 


-0.153 


Sen 


361.3840 


0.022 


-0.024 


Co i 


341.2338 


0.514 


-0.064 


Co i 


384.5468 


0.923 


-0.133 


















Co i 


384.5480 


0.923 


-0.227 


Sen 


424.6832 


0.315 


-0.096 


Co i 


344.9178 


0.581 


-0.174 


Co i 


384.5474 


0.923 


-0.143 


Sen 


424.6836 


0.315 


-0.080 


Co i 


344.9172 


0.581 


-0.149 


Co i 


384.5465 


0.923 


-0.116 


Sen 


424.6836 


0.315 


-0.039 


Co i 


344.9183 


0.581 


-0.149 


Co i 


384.5481 


0.923 


-0.227 


Sen 


424.6839 


0.315 


-0.081 


Co i 


344.9176 


0.581 


-0.339 


Co i 


384.5473 


0.923 


-0.138 


Sen 


424.6839 


0.315 


-0.078 


Co i 


344.9167 


0.581 


-0.129 


Co i 


384.5462 


0.923 


-0.101 


Sen 


424.6841 


0.315 


-0.098 


Co i 


344.9183 


0.581 


-0.129 


Co i 


384.5482 


0.923 


-0.237 


Sen 


424.6842 


0.315 


-0.148 


Co i 


344.9173 


0.581 


-0.215 


Co i 


384.5471 


0.923 


-0.139 


Sen 


424.6843 


0.315 


-0.291 


Co i 


344.9161 


0.581 


-0.122 


Co i 


384.5458 


0.923 


-0.088 


Sen 


424.6843 


0.315 


-0.096 


Co i 


344.9182 


0.581 


-0.122 


Co i 


384.5483 


0.923 


-0.257 


Sen 


424.6844 


0.315 


-0.116 


Co i 


344.9170 


0.581 


-0.143 


Co i 


384.5470 


0.923 


-0.146 


Sen 


424.6845 


0.315 


-0.080 


Co i 


344.9154 


0.581 


-0.125 


Co i 


384.5455 


0.923 


-0.077 


Sen 


424.6846 


0.315 


-0.098 


Co i 


344.9181 


0.581 


-0.125 


Co i 


384.5483 


0.923 


-0.297 


Sc ii 


424.6846 


0.315 


-0.081 


Co i 


344.9165 


0.581 


-0.103 


Co i 


384.5468 


0.923 


-0.167 










Co i 


344.9146 


0.581 


-0.144 


Co i 


384.5450 


0.923 


-0.066 


Co i 


340.5074 


0.432 


-0.139 


Co i 


344.9178 


0.581 


-0.144 










Co i 


340.5072 


0.432 


-0.159 


Co i 


344.9159 


0.581 


-0.074 


Co i 


399.5299 


0.923 


-0.203 



Co i 


340.5081 


0.432 


-0.159 










Co i 


399.5304 


0.923 


-0.203 


Co i 


340.5079 


0.432 


-0.138 


Co i 


345.3472 


0.432 


-0.105 


Co i 


399.5302 


0.923 


-0.305 


Co i 


340.5075 


0.432 


-0.137 


Co i 


345.3479 


0.432 


-0.150 


Co i 


399.5297 


0.923 


-0.175 


Co i 


340.5088 


0.432 


-0.137 


Co i 


345.3475 


0.432 


-0.093 


Co i 


399.5306 


0.923 


-0.175 


Co i 


340.5085 


0.432 


-0.124 


Co i 


345.3489 


0.432 


-0.250 


Co i 


399.5302 


0.923 


-0.235 


Co i 


340.5081 


0.432 


-0.126 


Co i 


345.3485 


0.432 


-0.129 


Co i 


399.5296 


0.923 


-0.161 


Co i 


340.5098 


0.432 


-0.126 


Co i 


345.3480 


0.432 


-0.081 


Co i 


399.5308 


0.923 


-0.161 


Co i 


340.5093 


0.432 


-0.107 


Co i 


345.3498 


0.432 


-0.229 


Co i 


399.5302 


0.923 


-0.190 


Co i 


340.5088 


0.432 


-0.122 


Co i 


345.3493 


0.432 


-0.119 


Co i 


399.5294 


0.923 


-0.154 


Co i 


340.5110 


0.432 


-0.122 


Co i 


345.3488 


0.432 


-0.070 


Co i 


399.5310 


0.923 


-0.154 


Co i 


340.5104 


0.432 


-0.090 


Co i 


345.3510 


0.432 


-0.226 


Co i 


399.5302 


0.923 


-0.157 


Co i 


340.5097 


0.432 


-0.123 


Co i 


345.3504 


0.432 


-0.115 


Co i 


399.5292 


0.923 


-0.153 


Co i 


340.5123 


0.432 


-0.123 


Co i 


345.3497 


0.432 


-0.059 


Co i 


399.5313 


0.923 


-0.153 


Co i 


340.5116 


0.432 


-0.074 


Co i 


345.3524 


0.432 


-0.234 


Co i 


399.5302 


0.923 


-0.131 


Co i 


340.5108 


0.432 


-0.131 


Co i 


345.3517 


0.432 


-0.117 


Co i 


399.5290 


0.923 


-0.160 


Co i 


340.5139 


0.432 


-0.131 


Co i 


345.3508 


0.432 


-0.050 


Co i 


399.5315 


0.923 


-0.160 


Co i 


340.5131 


0.432 


-0.060 


Co i 


345.3540 


0.432 


-0.253 


Co i 


399.5302 


0.923 


-0.109 


Co i 


340.5121 


0.432 


-0.152 


Co i 


345.3531 


0.432 


-0.125 


Co i 


399.5288 


0.923 


-0.180 


Co i 


340.5156 


0.432 


-0.152 


Co i 


345.3522 


0.432 


-0.041 


Co i 


399.5317 


0.923 


-0.180 


Co i 


340.5147 


0.432 


-0.046 


Co i 


345.3558 


0.432 


-0.292 


Co i 


399.5302 


0.923 


-0.090 










Co i 


345.3548 


0.432 


-0.147 










Co i 


340.9155 


0.514 


-0.204 


Co i 


345.3537 


0.432 


-0.032 


Co i 


412.1330 


0.923 


-0.213 


Co i 


340.9159 


0.514 


-0.204 










Co i 


412.1332 


0.923 


-0.206 


Co i 


340.9157 


0.514 


-0.306 


Co i 


349.5710 


0.629 


-0.147 


Co i 


412.1326 


0.923 


-0.186 


Co i 


340.9154 


0.514 


-0.176 


Co i 


349.5703 


0.629 


-0.147 


Co i 


412.1337 


0.923 


-0.276 


Co i 


340.9162 


0.514 


-0.176 


Co i 


349.5693 


0.629 


-0.177 


Co i 


412.1331 


0.923 


-0.186 


Co i 


340.9159 


0.514 


-0.236 


Co i 


349.5707 


0.629 


-0.165 


Co i 


412.1323 


0.923 


-0.166 


Co i 


340.9155 


0.514 


-0.162 


Co i 


349.5697 


0.629 


-0.131 


Co i 


412.1338 


0.923 


-0.260 


Co i 


340.9167 


0.514 


-0.162 


Co i 


349.5684 


0.629 


-0.134 


Co i 


412.1329 


0.923 


-0.176 


Co i 


340.9163 


0.514 


-0.191 


Co i 


349.5702 


0.629 


-0.190 


Co i 


412.1318 


0.923 


-0.149 


Co i 


340.9157 


0.514 


-0.155 


Co i 


349.5689 


0.629 


-0.130 


Co i 


412.1338 


0.923 


-0.260 


Co i 


340.9173 


0.514 


-0.155 


Co i 


349.5673 


0.629 


-0.105 


Co i 


412.1327 


0.923 


-0.171 


Co i 


340.9167 


0.514 


-0.158 


Co i 


349.5695 


0.629 


-0.233 


Co i 


412.1313 


0.923 


-0.134 


Co i 


340.9160 


0.514 


-0.154 


Co i 


349.5679 


0.629 


-0.146 


Co i 


412.1338 


0.923 


-0.270 


Co i 


340.9180 


0.514 


-0.154 


Co i 


349.5659 


0.629 


-0.084 


Co i 


412.1324 


0.923 


-0.172 


Co i 


340.9173 


0.514 


-0.132 










Co i 


412.1308 


0.923 


-0.121 


Co i 


340.9164 


0.514 


-0.161 


Co i 


350.2243 


0.432 


-0.174 


Co i 


412.1338 


0.923 


-0.290 


Co i 


340.9188 


0.514 


-0.161 


Co i 


350.2242 


0.432 


-0.167 


Co i 


412.1321 


0.923 


-0.179 


Co i 


340.9179 


0.514 


-0.110 


Co i 


350.2238 


0.432 


-0.237 


Co i 


412.1301 


0.923 


-0.110 


Co i 


340.9169 


0.514 


-0.181 


Co i 


350.2248 


0.432 


-0.147 


Co i 


412.1336 


0.923 


-0.330 


Co i 


340.9197 


0.514 


-0.181 


Co i 


350.2245 


0.432 


-0.147 


Co i 


412.1316 


0.923 


-0.200 


Co i 


340.9187 


0.514 


-0.091 


Co i 


350.2239 


0.432 


-0.221 


Co i 


412.1294 


0.923 


-0.099 










Co i 


350.2255 


0.432 


-0.127 










Co i 


341.2325 


0.514 


-0.178 


Co i 


350.2249 


0.432 


-0.137 











